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PREFACE 


A  geophysical  Investigation  of  foundation  conditions  of  the  right  abutment 
area  of  Mill  Creek  Dam,  Walla  Wallat  Washington »  was  authorized  by  the  US  Army 
Engineer  District,  Walla  Walla  (CENPW) ,  under  lAO  No.  E86890074,  dated  4  April 
1989.  The  work  was  performed  during  the  period  4  April- 30  September  1989. 

Program  planning,  field  work,  data  processing,  and  Interpretation  was 
accomplished  by  a  team  consisting  of  Messrs.  Donald  E.  Yule,  Michael  K.  Sharp, 
Jose  L.  Llopls  and  Dr.  Dwaln  K.  Butler,  Engineering  Geophysics  Branch  (EGB) , 
Earthquake  Engineering  and  Geosciences  Division  (EEGD) ,  Geotechnical  Laboratory 
(GL) ,  US  Army  Engineer  Waterways  Experiment  Station  (WES).  Personnel  of  CENPW 
performed  location  and  elevation  surveys  for  the  geophysical  survey  lines. 

This  work  was  performed  under  the  general  supervision  of  Mr.  Joseph  R.  Curro, 
Chief,  EGB,  Dr.  Arley  G.  Franklin,  Chief,  EEGD,  and  Dr.  William  F.  Marcuson, 
Chief,  GL.  This  report  was  prepared  by  Dr.  Butler. 

COL  Larry  B.  Fulton,  EN,  was  Commander  und  Director  of  WES  during  the 
publication  of  this  report.  Dr.  Robert  W.  Whalln  was  Technical  Director. 
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GEOPHYSICAL  ASSESSMENT  OF  FOUNDATION  CONDITIONS 
RIGHT  ABUTMENT,  MILL  CREEK  DAM 


PART  I  :  INTRODUCTION 


Background 

1.  Mill  Creek  Dam,  near  Walla  Walla,  Washington,  has  ex-* 
perienced  anomalous  seepage  since  its  first  filling  in  1941. 
Various  attempts  to  abate  and  control  the  seepage,  including  con¬ 
struction  of  a  concrete  cutoff  wall  (completed  in  1983),  have 
been  unsuccessful.  As  part  of  a  test  pool  raise  in  1984,  WES  in¬ 
stalled  a  self  potential  array  to  monitor  seepage  conditions 
during  the  test.  The  results  of  the  geophysical  surveys  and 
other  geotechnical  investigations  indicate  that  the  majority  of 
the  anomalous  seepage  is  presently  occurring  through  the  right 
abutment  of  the  dam,  around  the  end  of  the  cutoff  wall  (Butler, 
Wahl  and  Sharp,  1984) .  The  geophysical  surveys  also  indicated  a 
possible  minor  seepage  path  under  the  embankment  near  the  left 
abutment,  where  a  deficiency  may  exist  in  the  cutoff  wall  (the 
cutoff  wall  may  not  extend  deep  enough).  Overall,  the  seepage 
has  been  reduced  by  30  percent,  from  33  ftvs  (0.93  m^/s)  to  22 
ft^/s  (0.62  icr/s)  ,  since  construction  of  the  cutoff  wall,  and 
downstream  flooded  (saturated)  farmland  acreage  has  been  reduced 
by  56  percent.  However,  seepage  pressures  in  a  conglomerate  for¬ 
mation  in  the  right  abutment  area  have  increased  dramatically 
(Walla  Walla  District,  1987) . 

2.  The  dam  foundation  and  abutments  are  silt  (the  Palouse 
silt  or  loess) ,  which  varies  from  50  to  100  ft  thick  beneath  the 
upstream  toe  of  the  dam.  Below  the  silt  are,  successively,  the 
conglomerate  and  a  basalt.  The  conglomerate  (alluvial)  consists 
of  coarse  sand,  gravel,  cobbles,  and  zones  of  silt;  while  the 
conglomerate  in  the  vicinity  of  the  dam  has  no  formal  name,  it  is 
correlative  to  the  Ringold  Formation.  The  permeability  of  the 
conglomerate  varies  considerably  both  laterally  and  with  depth  in 
the  formation.  During  the  initial  and  subsequent  reservoir  fill¬ 
ing  events,  sinkholes  and  cracks  appeared  in  the  reservoir  floor 
and  right  abutment  area.  Also,  significant  quantities  of  embank¬ 
ment  and/or  foundation  materials  were  piped  through  the  dam^s  in¬ 
terior  drainage  system.  The  project  history  coupled  with  the  al¬ 
tered  seepage  conditions  caused  by  the  presence  of  the  cutoff 
wall  lead  to  concerns  about  the  safety  of  the  dam  during  subse¬ 
quent  reservoir  filling.  A  failure  scenario  for  the  dam  is  based 
on  a  model  which  includes  piping  of  the  silt 
foundation/embankment  materials  through  the  conglomerate  leading 
to  subsidence  of  overlying  materials.  Formation  and  growth  of 
voids  and  subsidence  in  the  foundation/ embankment  could  threaten 
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the  Integrity  of  the  dam.  The  basalt,  the  Frenchman's  Springs 
Member  of  the  Wanatum  Basalt  (Columbia  River  Basalt  Group) ,  is 
highly  fractured  in  places  and  has  portions  which  are  vesicular. 

3.  Walla  Walla  District  has  proposed  an  Impermeable 
membrane  liner  for  the  reservoir.  Prior  to  approval  of  the 
design  concept,  additional  foundation  investigations  including 
surface  geophysical  surveys  were  planned.  WES  was  asked  to  plan 
a  geophysical  field  program  which  concentrated  on  assessing  foun¬ 
dation  conditions  in  the  right  abutment  area. 


Scope 

4.  WES  designed  a  geophysical  field  survey  program  to 
(1)  assess  general/overall  subsurface  conditions  in  the  right 
abutment  area,  such  as  depth  to  formation  interfaces  and  lateral 
condition  changes,  and  (2)  detect  localized  anomalous  conditions, 
such  as  cavities,  subsidence  features,  saturated  zones,  and  any 
anomalous  low  density  zones.  Part  II  of  this  report  describes 
the  field  survey  program  and  rationale  and  presents  a  brief  sum¬ 
mary  of  each  geophysical  method  employed  at  the  site.  Part  III 
presents  the  results  and  interpretation  of  the  geophysical 
program,  and  a  summary  and  recommendations  are  given  in  Part  IV. 
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PART  II  :  GEOPHYSICAL  SURVEY  PROGRAM 


Program  Design 


Site  Speeifio  Details 

5.  The  stratigraphy  of  the  foundation  and  right  abutment 
is  known  to  the  extent  illustrated  in  Figure  1  (from  a  Draft 
Design  Memorandum,  Walla  Walla  District,  1987),  nowever  the  dis¬ 
tribution  of  physical  properties  of  the  silt  and  conglomerate  is 
uncertain.  The  permeability  of  the  conglomerate,  which  consists 
of  coarse  sand,  gravel,  cobbles,  and  zones  of  silt,  for  example, 
is  known  to  vary  considerably.  The  variable  permeability  of  the 
conglomerate  is  directly  related  to  the  degree  of  cementation. 
Also,  the  nature  of  anomalous  conditions,  if  they  exist  at  all, 
implied  by  the  failure  scenario  is  completely  uncertain.  Piping 
of  the  silt  into  the.  underlying  conglomerate  could  produce  voids 
in  the  silt  formation,  but  more  likely  the  piping  would  result  in 
a  progressive  collapse  feature  which  would  migrate  upward  and 
laterally  away  from  the  location  where  the  silt  enters  the  con¬ 
glomerate.  Within  the  silt  then,  the  effect  would  be  to  produce 
anomalously  low  density  zones.  Within  the  conglomerate,  the 
anomalous  seepage  would  carry  downstream  not  only  the  material 
removed  from  the  silt  above  but  material  removed  from  within  the 
conglomerate  itself,  such  as  fines  from  within  the  matrix  of 
uncemented  zones  of  conglomerate  or  from  silt  pockets  in  the 
conglomerate.  The  net  effect  would  be  to  produce  zones  of 
anomalously  low  density  and/or  higher  water  content  within  the 
conglomerate . 

6.  The  most  obvious  physical  effect  of  the  piping/collapse 
model  will  be  the  creation  of  low  density  zones  in  the  silt  end 
conglomerate.  Electrical  resistivity  and  conductiv  ity 
(resistivity  -  1  /  conductivity)  depend  intimately  on  porosity 
and  water  content  (saturation) .  The  piping  should  produce  zones 
of  increased  porosity  and,  where  the  zones  are  below  the  water 
table,  increased  water  content.  Also,  the  zones  of  increased 
porosity  may  retain  water  from  rainfall  infiltration  and 
anomalous  seepage.  The  higher  porosity/ water  content  zones  will 
produce  low  resistivity/high  conductivity  anomalous  zones.  Where 
the  anomalous  zones  in  the  silt  are  saturated,  high  seismic 
velocity  zones, will  be  created. 

7.  It  was  difficult  in  the  geophysical  program  planning 
phase  to  optimize  a  program  for  the  site,  since  nothing  was  known 
about  such  geophysical  parameters  as  density,  electrical  resis¬ 
tivity,  and  seismic  velocity  of  the  silt,  conglomerate  and 
basalt.  Also,  as  discussed  above,  the  details  of  anomalous  con¬ 
ditions  associated  with  the  anomalous  seepage  are  unknown.  Thus 
the  geophysical  program  had  to  be  planned  to  include  a  range  of 
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in  situ  condition  possibilities.  The  strategy  was  to  apply  com¬ 
plementary  geophysical  methods  and  to  achieve  high  resolution 
subsurface  coverage  both  vertically  and  laterally. 


Complementary  Geophysical  Program 

8.  The  geophysical  program  includes  the  following  methods: 
seismic  reflection  (SR) ,  ground  penetrating  radar  (GPR) , 
electromagnetic  conductivity  (EM) ,  electrical  resistivity  (ER)  , 
and  microgravity  (MG) .  If  site  conditions  are  such  as  to  allow 
sufficient  depth  of  investigation  (penetration) ,  the  GPR  method 
holds  the  potential  for  continuous  lateral  coverage  and  the 
highest  resolution  of  any  of  the  geophysical  methods.  Depth  of 
investigation  for  GPR  varies  greatly  from  site  to  site,  from  as 
much  as  50  ft  to  less  than  3  ft  depending  on  the  ground  conduc¬ 
tivity  (penetration  decreases  as  conductivity  increases) .  Clays, 
high  water  content  and  shallow  water  tables  are  site  factors 
which  can  limit  GPR  penetration.  Localized  features  as  small  as 
2-3  ft  effective  diameter  can  generally  be  detected  within  the 
effective  depth  of  investigation,  if  the  features  are  conductive 
or  highly  resistive  relative  to  the  surrounding  material .  GPR  is 
also  highly  effective  for  profiling  interface  depth  variations. 
Seismic  reflection  complements  GPR  in  providing  relatively  high 
resolution  laterally  and  vertically,  however  SR  has  consistently 
greater  depth  of  investigation  capability.  SR  has  the  potential 
for  detecting  localized  features  ar  small  as  10-15  ft  in  effec¬ 
tive  diameter.  Objectives  of  the  SR  and  GPR  surveys  are  to  map 
interfaces  between  the  silt/conglomerate  and  conglomerate/basalt 
and  to  detect  localized  anomalous  conditions  such  as  voids  or 
collapse  structures. 

9.  EM  and  ER  are  complementary  in  the  sense  that  the 
methods  respond  with  different  sensitivities  to  high  and  low 
resistd vity/conductivity  anomalies.  EM  can  be  used  to  profile 
long  dijtances  rapidly  and  is  sensitive  to  high  conductivity  lo¬ 
calized  features.  ER  is  well  suited  to  both  vertical  property 
variation  detection  and  horizontal  profiling,  and  is  ap¬ 
proximately  equally  sensitive  to  high  and  low  resistivity 
anomalies.  ER  surveying  is  more  labor  intensive  than  EM  survey¬ 
ing.  ER  and  EM  survey  objectives  are  to  map  variations  in 
material  properties  such  as  porosity  and  water  content  and  to 
detect  localized  anomalies  such  as  voids  and  collapse  features. 

10.  Microgravity  surveying  gives  a  positive  indication  of 
low  and  high  density  anomalous  conditions  in  the  subsurface. 
Coupled  with  the  EM  and  ER  results,  the  results  of  microgravity 
surveys  allow  a  higher  probability  anomaly  identification:  such 
as,  air-filled  versus  saturated  higher  porosity  zone  (cavity  or 


collapse  feature)  ;  cemented,  low  porosity/permeability  con** 
glomerate  versus  unuemented,  high  porosity/permeabllity  con-* 
glomerate. 

Program  Layout 

11.  Figure  2  is  a  site  location  map  of  the  right  abutment 
showing  geophysical  survey  lines.  The  survey  lines  were  estab¬ 
lished  to  provide  continuous  coverage  along  the  length  of  the  em¬ 
bankment  and  approximately  perpendicular  to  the  axis  of  the  em¬ 
bankment  in  the  right  abutment  area.  Survey  lines  C-F  begin  at  a 
location  approximately  240  ft  from  the  end  of  the  cutoff  wall 
near  station  32+20.  Walla  Walla  District  survey  personnel  lo¬ 
cated  the  lines  and  determined  elevations  and  coordinates  of  sur¬ 
vey  points.  For  lines  C-F,  elevations  were  determined  every 
10  ft  and  coordinates  every  20  ft.  Stations  are  numbered  con¬ 
secutively  along  each  survey  line,  beginning  from  the  end  nearest 
the  cutoff  wall.  Odd-numbered  stations  are  indicated  by  an  off¬ 
set  numbered  stake  and  a  wooden  hub  or  steel  spike  driven  flush 
with  the  ground  L.urface.  Even-numbered  stations  are  indicated  by 
a  paint  spot.  Along  the  crest,  line  B  begins  at  borehole  RDH-6, 
and  stations  are  located  every  20  ft  along  the  center  of  the  ac¬ 
cess  road. 

12.  In  addition  to  the  profile  lines  in  Figure  2,  a  survey 
grid  was  established  around  a  sinkhole  as  illustrated  in 
Figure  3 .  The  sinkhole  is  located  in  the  floor  of  the  reservoir 
just  off  the  Figure  2  site  map,  200  ft  NE  of  the  end  of  line  D 
(D59)  and  250  ft  nearly  due  east  the  end  of  line  F  (P65) .  The 
sinkhole  survey  grid  was  established  and  elevations  determined  by 
WES  personnel.  The  sinkhole  is  approximately  10  ft  in  diameter 
and  3  ft  deep.  Apparently  the  present  sinkhole  is  in  the 
vicinity  of  a  previous  sinkhole  which  was  filled. 

13 .  The  complete  geophysical  program  plan  is  summarized  in 
Table  1,  which  indicates  the  survey  line  length,  measurement  sta¬ 
tions,  and  the  geophysical  methods  applied  along  each  line. 


Geophysical  Methods  Concents  and  Field  Procedures 

14.  No  attempt  will  be  made  to  thoroughly  review  the 
physical  concepts  of  the  geophysical  methods.  For  each  method, 
a  very  brief  concept  statement  is  given  and  the  specifics  of  the 
field  procedure  for  this  site  presented.  For  each  method, 
reference  will  be  made  to  specific  readily  available  WES  techni¬ 
cal  reports  or  other  appropriate  sources.  General  references  for 
geophysical  methodology  are  Department  of  the  Army  (1979)  and 
Telford  et  al  (1976). 
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N-S  SURVEY  LINES 


Figure  3.  Survey  grid  around  sinkhole.  Elevation  contour  interval 


Table  1 


Geophysical  Program 


Survey  Line 

Line/Area  Length  f Stations^  Geophysical  Methods- 


B 

600  ft 
(B1-B31) 

MG 

C 

740  ft 
(C1-C75) 

SR,  GPR,  ER,  EM, 

D 

580  ft 
(D1-D59) 

GPR,  SR,  MG 

E 

300  ft 
(E1-E31) 

GPR,  EM,  MG 

F 

640  ft 
(F1-F65) 

GPR,  EM,  MG 

Sinkhole 

100  X  100  ft 
survey  grid 

GPR,  EM,  MG 

*MG — Microgravity?  SR — seismic  reflection;  GPR — ground 
penetrating  radar;  lR--electrical  resistivity; 

EM — electromagnetic  conductivity. 
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15.  The  seiemle  reflection  method  has  been  a  cornerstone 
of.  petroleum  exploration  for  decades,  and  recent  advances  in 
seismic  instrumentation  and  microcomputers  make  it  feasible  to 
apply  the  same  methodology  to  geotechnical  problems  (shallow 
depths  of  investigation,  from  as  shallow  as  a  few  feet  to  several 
hundred  feet) .  Seismic  reflection  involves  near  vertical 
propagation  paths  in  the  subsurface,  contrasted  to  the  seismic 
refraction  method  which  utilizes  horizontal  travel  paths,  and 
seeks  to  identify  reflections  from  subsurface  interfaces  or  lo¬ 
calized  anomalous  features.  The  practical  implication  of  the 
near  vertical  propagation  paths  is  that  the  distances  from  the 
seismic  source  to  the  farthest  geophone  is  much  shorter  than  for 
seismic  refraction  for  the  same  depth  of  investigation.  Thus 
seismic  reflection  can  be  effectively  applied  at  sites  with 
limited  horizontal  access,  such  as  dams  sites.  Dobecki  et  al 
(1989)  discuss  the  seismic  reflection  method  and  its  application 
at  a  Corps'  dam  site. 

16.  Two  different  24-channel  seismographs  were  tested. 
Three  seismic  sources,  a  20-lb  sledge  hammer,  a  downhole  30-06 
rifle,  and  a  downhole  .50  caliber  rifle,  were  tested  at  several 
different  filter  settings  of  the  seismographs.  The  field 
parameter  selection  array  (geophone  layout)  consisted  of  24 
takeouts  at  2-ft  spacing,  where  each  takeout  was  connected  to 
three  closely  spaced  40-Hz  geophones  in  series.  Field  parameter 
selection  tests  were  conducted  at  three  locations,  two  locations 
on  line  C  and  one  location  on  line  D.  At  each  location  the  in¬ 
dividual  sources  were  tested  at  various  filter  settings  and  at 
offset  distances  ranging  from  2  to  96  ft  from  the  closest 
geophone.  This  parameter  selection  test  procedure  is  termed  a 
walkaway  test.  For  each  source  and  filter  setting  tested,  the 
walkaway  test  resulted  in  records  covering  either  96  or  144  ft 
laterally.  The  work  was  performed  by  Great  Plains  Geophysical, 
Lawrence ,  Kansas . 

Ground  Penetrating  Radar 

17.  The  GPR  method  produces  a  continuous  profile  of 
electromagnetic  (radar)  reflections  from  subsurface  features. 
This  is  accomplished  by  towing  a  radar  antenna  at  a  walking  pace 
along  the  desired  survey  line.  The  antenna  is  connected  to  a 
controller/processor  and  graphic  recorder  mounted  on  a  cart  or 
field  vehicle.  During  the  traverse  along  the  survey  line,  the 
antenna  alternately  transmits  and  receives  radar  pulses  (120  MHz 
center  frequency) .  The  received  signals  represent  reflections 
(near  vertical  propagation  paths)  from  subsurface  features.  The 
received  signals  are  displayed  on  the  graphic  recorder  as  a 
record  of  signal  amplitude  versus  time  ("radargram”) -  Figure  4 
illustrates  the  GPR  technique  (Sylvester,  1989).  The  pulse  is 
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CABLE  TO  OSPLAY  RECOBOER 


b.  Graphic  radar  record  ("radargram")  for  geologic  section  in  a 


Figure  4.  Illustration  of  GPR  technique  and  graphic  radar  record. 
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rep<iated  and  the  antenna  is  towed  at  a  rate  which  gives  3-4  ver¬ 
tical  records  (soundings)  every  foot  along  the  survey  line. 
Along  lines  C-F,  records  were  obtained  using  100  and  300 
nanoseconds  (ns)  total  record  tine,  l.e.,  two  traverses  were  per¬ 
formed  along  each  line.  Twelve  survey  lines  were  traversed  in 
the  area  of  the  sinkhole)  record  time  was  150  ns,  except  for  line 
1  at  200  ns.  As  the  antenna  Is  towed,  manually  triggered 
horizontal  distance  Indicator  marks  are  produced  on  the  record  as 
each  survey  stake  Is  passed.  Ballard  (1983)  discusses  the  con¬ 
cepts  of  the  GPR  method  and  presents  examples  of  6PR  records  over 
subsurface  features  such  as  cavities.  The  GPR  work  was  performed 
by  Williamson  and  Associates,  Seattle,  Washington. 


Electromagnetic  Conductivity 

18 .  The  EM  method  Involves  the  determination  of  an  ap¬ 
parent  ground  conductivity  from  measurements  with  a  system  that 
consists  of  an  electromagnetic  loop  transmitter  and  a  loop 
receiver.  The  transmitter  emits  a  constant  frequency  signal  (the 
primary  signal)  at  400,  1600  or  6400  Hz.  The  receiver  detects  a 
signal  that  consists  of  the  primary  signal  plus  a  secondary  sig¬ 
nal.  The  secondary  signal  arises  from  Interactions  with  subsur¬ 
face  materials  and  features  through  the  process  of  electromag¬ 
netic  induction.  A  ratio  of  the  secondary  signal  to  the  primary 
signal  Is  then  formed  by  the  receiver  electronics.  The  ratio  Is 
proportional  to  the  apparent  ground  conductivity,  and  the 
receiver  displays  a  reading  In  conductivity  units  (millimhos/m  » 
10“^/ohm-m,  where  ohm-m  are  the  units  of  resistivity) .  If  the 
subsurface  were  a  uniform  material  to  a  great  depth,  the  apparent 
ground  conductivity  would  equal  the  true  conductivity  of  the 
material.  In  general  however,  the  conductivity  reading  Is  a 
weighted  average  of  all  material  within  the  volume  of  investiga¬ 
tion.  The  volume  of  investigation  and  hence  the  depth  of  Inves¬ 
tigation  depends  primarily  on  the  transmitter  and  receiver  loop 
separation  and  the  transmitter  frequency.  The  basic  concept  of 
EM  surveying  is  illustrated  in  Figure  5. 

19.  Loop  separations  of  10  m,  20  m,  and  40  m  are  used 
respectively  for  the  three  transmitter  frequencies  listed  above. 
For  each  loop  separation,  measurements  can  be  made  for  two  loop 
orientations;  horizontal  loops  (vertical  dipoles,  VD) ,  and  verti¬ 
cal  loops  (horizontal  dipoles,  HD) .  Thus  at  each  location  six 
different  measurements  can  be  made.  The  approximate  depths  of 
investigation  for  the  six  measurements  are  listed  below: 
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m  (ft). 

m 

m 

10 

7.5 

15 

(33) 

(25) 

(50) 

20 

15 

30 

(66) 

(50) 

(100) 

40 

30 

60 

(130) 

(100) 

(200) 

within  the  approximate  depths  of  investigation  listed  above,  the 
HD  measurements  are  the  most  sensitive  to  shallow  material.  All 
six  measurements  were  made  at  lO  ft  intervals  along  lines  C,E  and 
F.  The  transmitter  was  placed  at  the  survey  line  station  for 
each  set  of  measurements.  Thus,  following  the  convention  of  as¬ 
signing  the  measurement  to  the  midpoint  of  the  loop  separation, 
measurements  with  the  three  loop  separations  will  be  displaced 
different  amounts  from  each  survey  line  station  in  data  plots. 
Butler  (1986)  discusses  the  EM  method  and  presents  examples. 

Electrical  Resistivity 

20.  ER  methods  are  used  generally  to  determine  vertical 
(vertical  sounding)  or  lateral  (horizontal  profiling)  profiles  of 
electrical  resistivity  of  the  subsurface.  The  method  used  for 
the  present  investigation  utilizes  a  pole-dipole  electrode  array 
and  is  actually  a  combined  vertical  sounding/horizontal  profiling 
method.  Field  procedure  and  interpretation  method  are  a 
modification  by  Bates  (1973)  of  a  methodology  first  proposed  by 
Bristow  (1966),  Details  of  field  procedures,  data  processing 
methods,  and  interpretation  concepts  are  given  by  Butler  et  al 
(1982) .  The  major  disadvantages  of  the  pole-dipole  method  of 
Bristow  and  Bates  are  the  very  labor  intensive  data  acquisition 
procedure  and  the  subjective  interpretation  procedure.  The  major 
advantages  are  that  the  interpretation  procedure  locates 
anomalies  in  cross  section  and  gives  an  indication  of  size  and 
shape.  Numerous  case  studies  indicate  successful  application  of 
the  method  to  detection  and  delineation  of  subsurface  cavities. 

21.  A  pole-dipole  resistivity  survey  was  conducted  along 
line  C.  Beginning  at  station  Cl,  current  electrode  stations  were 
established  every  50  ft  along  the  line.  A  second  current 
electrode  was  placed  approximately  1000  ft  distant  which  remained 
fixed  during  the  survey.  For  each  current  electrode  station 
along  the  line,  two  vertical  soundings  are  performed  by  moving  a 
potential  electrode  pair,  separated  by  10  ft,  along  the  profile 
line  from  10  ft  to  200  ft  on  each  side  of  the  current  electrode. 
The  resulting  data  set  is  extensive  and  has  considerable  redun- 
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dancy.  Tha  vertical  soundings  are  examined  and  vertical  varia-^ 
tlon  trends  are  determined.  For  each  sounding,  high  and  low 
resistivity  anomalies  relative  to  the  trend  are  Identified. 
Using  the  procedure  discussed  in  Butler  et  al  (1982)  the 
anomalies  are  located  on  a  cross  section  beneath  line  G. 
Anomalies  are  defined  in  terms  of  profile  location,  depth  and  ap¬ 
proximate  shape  and  size.  The  graphical  interpretation  technique 
is  Illustrated  in  Figure  6a,  and  an  example  of  the  sounding 
curves  and  Interpreted  cross-section  for  a  survey  across  an  air- 
filled  cavity  system  is  shown  in  Figure  6b, c  (Butler  et  al  1982) . 

Mlcroaravltv 

22.  The  microgravity  method  involves  relative  measurements 
of  the  acceleration  or  force  of  gravity  on  the  surface  with  high 
accuracy  and  precision.  A  gravity  meter  is  used  with  an  instru¬ 
ment  sensitivity  of  1  microgal,  where  1  mlcrogal  ==  (10"^  x  g)  and 
g  is  the  normal  earth's  gravitational  acceleration  of  980  cm/s^ 
(32  tt/wr) .  After  applying  various  corrections,  the  gravity  data 
are  presented  in  profile  plots  or  contour  plots.  If  the  subsur¬ 
face  is  uniformly,  horizontally  layered  beneath  the  survey  site, 
for  example,  the  measured  gravity  field  would  be  a  constant  or 
perhaps  be  slowly  varying  (termed  the  regional  gravity  field  com¬ 
ponent)  over  the  site  due  to  deep  geologic  conditions.  The 
presence  of  any  lateral  density  changes  at  shallow  depths  beneath 
the  site  will  produce  gravity  anomalies  which  will  perturb  the 
constant  or  slowly  varying  regional  gravity  field.  For  example, 
a  2  m  (6.6  ft)  diameter  spherical-shaped,  air-filled  cavity  in 
soil  can  be  detected  to  a  depth  of  about  4  m  (13  ft),  whereas  a 
2  m  diameter  horizontal,  cylindrical-shaped,  air-filled  cavity 
can  be  detected  to  a  depth  of  about  9m  (30  ft) .  Figure  7  il¬ 
lustrates  the  concept  of  a  gravity  anomaly  produced  by  a  subsur¬ 
face  density  anomaly  and  shows  gravity  anomalies  over  spherical 
and  cylindrical  features  for  three  depths.  Note  that  the  anomaly 
magnitude  decreases  as  depth  increases,  while  the  spatial  extent 
of  the  anomaly  increases  as  depth  increases. 

23.  Butler  (1980)  describes  the  field  procedures  and  data 
corrections  required  for  microgravity  surveying.  The  data  must 
be  acquired  using  a  very  careful  and  consistent  field  procedure 
which  involves  an  exacting  measurement  sequence,  repeated  obser¬ 
vations  at  20  percent  or  more  of  the  stations,  and  repeated  oc¬ 
cupations  of  base  stations  at  30-45  minute  intervals.  Data  cor¬ 
rections  Include: 

a.  Time  variations  of  measured  gravity  which  consist  of 
Instrument  drift  and  earth  tide  variations;  base  station 
reoccupations  are  utilized  for  this  correction 
procedure ; 
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•  .  Graphical  intarpretatlon  method  lor  pola-ulpola  raaiativity  survaya. 
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0.  Comparison  of  resistivity  cross-section  with  geologic  cross-section. 

Figure  6.  Illustration  of  the  pole-dipole  resistivity  interpretation  procedure 
and  a  field  example. 
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MEASUREMENT  STATIONS 


GEOMETRY  OF  SPHERICAL  MODEL 


HORIZONTAL  CYLINDRICAL  MODEL  GEOMETRY 
MODEL  GEOMETRY 


a.  Concept  of  gravity  meaauiementt  acroas  an  anomaly  A  donaity  anomaly  /^produces 
a  gravity  anomaly  ySl  The  gravity  motei  measuraa  the  vertical  component  09a» 


A|,  iiOt~ 


b  Gravity  anomalies  aerosa  a  spherical  model  lor  three  depths  Zo. 
6t,.  *io« 


X.  M 

c  Gravity  anomalies  across  a  cylindrical  model  lor  three  depths  2o. 


Figure  7.  Illustration  of  gravity  anomalies  produced  by  simple  models 
of  subsurface  features. 
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b.  Latitude  variations  of  gravity,  l.e.,  correction  for 
distance  north/ south  of  a  base  station; 

c.  Elevation  variation  along  survey  lines; 

d.  Terrain  surrounding  the  survey  area. 

The  terrain  correction  for  a  survey  site  such  as  Mill  Creek  Dam  s 
quite  Involved,  requiring  consideration  of  terrain  variations  to 
a  distance  of  several  kilometers;  although  It  Is  only  terrain 
features  within  several  hundred  meters  of  a  measurement  station 
which  must  be  modeled  carefully.  Gravity  and  topographic  data 
from  the  survey  were  taken  to  the  Branch  of  Geophysics,  U.S. 
Geological  Survey,  Denver,  for  terrain  correction  calculations 
using  a  mainframe  computer  program.  The  remainder  of  the  gravity 
corrections  were  made  using  microcomputer  programs  described  by 
Butler  and  Yule  (1984) .  Gravity  data  which  has  been  corrected 
for  the  above  considerations  is  termed  the  Bouguer  gravity.  If 
the  constant  or  slowly  varying  regional  gravity  field  Is  sub¬ 
tracted  from  the  Bouguer  gravity,  the  remainder  Is  termed  the 
residual  gravity,  and  anomalies  are  termed  respectively  Bouguer 
and  residual  anomalies. 
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PAtlT  III  :  RESULTS  AND  INTEGRATED  ASSESSMENT 


Program  Results 


Seismic  Reflection 

24.  A  total  of  101  field  parameter  selection  tests  were 
conducted  at  the  three  locations  (source  locations  at  C13,  C41, 
and  D35) .  The  objectives  of  the  tests  were  to  select  source 
type,  geophone  type  (natural  frequency)  and  arrangement,  and 
seismograph  filter  settings  In  order  to  optimize  and  enhance  the 
detection  and  Identification  of  reflections.  Generally,  if  the 
shallow  SR  method  Is  suitable  for  the  geotechnical  objectives  at 
a  site,  reflection  events  can  be  identified  and  enhanced  during 
this  test  phase.  The  optimum  set  of  field  parameters  will  then 
be  used  for  the  remainder  of  the  reflection  survey  unless  condi¬ 
tions  change  significantly. 

25.  In  all  of  the  test  records,  only  one  event  is  noted 
which  is  considered  to  be  a  reflection.  The  reflection  event  is 
indicated  in  Figure  8  in  a  record  section  for  location  C13.  This 
event  has  a  normal  moveout  velocity  of  3730  ft/s  (1137  m/s)  and 
represents  reflection  from  an  interface  at  a  depth  of  170  ft. 
The  normal  moveout  velocity  represents  an  "average”  of  the  seis¬ 
mic  velocities  of  all  the  material  above  the  interface.  The 
reflection  is  probably  from  the  conglomerate-basalt  contact,  and 
the  depth  is  consistent  with  drilling  results  (Walla  Walla  Dis¬ 
trict  1987).  No  reflection  event  from  the  silt-conglomerate  con¬ 
tact  is  observed.  The  lack  of  identifiable  reflections  indicates 
that  the  acoustic  impedance  (product  of  bulk  density  and  seismic 
velocity)  contrast  across  the  silt-conglomerate  contact  must  be 
small.  This  lack  of  acoustic  impedance  contrast  likely  indicates 
that  the  conglomerate  is  uncemented  at  the  three  locations 
tested.  Also,  the  contact  may  be  gradational  and  irregular 
laterally;  both  conditions  would  tend  to  degrade  the  quality  of 
the  SR  records.  Figures  9  and  10  present  typical  records  from 
locations  C41  and  D35,  and  the  contractor's  report  is  included  as 
Appendix  A. 


26.  Refraction  events  are  present  on  all  records  at  the 
three  locations;  however,  due  to  the  short  seismic  survey  line 
lengths,  the  ii^ef factions  are  from  interfaces  (material  type  or 
condition  changes)  shallower  than  the  silt-conglomerate  contact. 
The  following  tabulation  summarizes  the  refraction  interpreta¬ 
tions  for  the  three  locations: 
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Figure  8.  Walkaway  seismic  test,  line  C,  location  C13. 
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Figure  9,  Walkaway  seismic  test,  line  C,  location  C41. 
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Figure  10.  Walkaway  seismic  test,  line  D,  location  D35. 
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Thus  the  compression  wave  (P-wave)  velocity  range  of  the  surface 
layer  (•'soil”)  is  650  -  1500  ft/s  (198  -  457  m/s)  and  the 

velocity  range  of  the  underlying  material  (silt/loess  with  vary¬ 
ing  water  content  and  calcite  content)  is  1400  -  4800  ft/s  (427  - 
1463  m/s)  .  The  4800  ft/s  velocity  at  the  first  site  (C13)  is 
clearly  anomalous,  and  likely  represents  a  refraction  off  the  top 
of  the  concrete  cutoff  wall.  The  maximum  interface  depth  at 
which  a  second  Interface  refraction  might  have  been  detected  by 
the  seismic  surveys  is  approximately  50  ft  (15  m) . 

Ground  Penetrating  Radar 

27.  Due  to  the  relatively  high  conductivity/low  resis¬ 
tivity  of  the  near  surface  materials  at  the  site,  the  maximum 
depth  from  which  radar  reflections  are  evident  is  20-30  ft.  The 
GPR  surveys  along  lines  C-P  indicate  predominately  uniform  condi¬ 
tions  and  thus  do  not  indicate  many  shallow  anomalies.  While  the 
GPR  results  are  not  useful  for  mapping  the  silt-conglomerate  in¬ 
terface  or  for  detecting  cavities  near  that  interface,  it  is  pos¬ 
sible  that  the  shallower  events  on  the  radar  records  could  be  in¬ 
dicative  of  deeper  anomalous  conditions.  The  GPR  surveys  in  the 
vicinity  of  the  sinkhole,  however,  contain  a  tremendous  amount  of 
information  indicating  subsurface  variation/anomalous  conditions 
in  the  upper  20-30  ft. 

28.  The  contractor's  report  is  included  as  Appendix  B  and 
contains  all  graphic  radar  records.  Only  representative  records 
will  be  presented  in  this  section.  Figure  11  is  the  record  for 
the  GPR  survey  along  line  E  and  is  typical  of  the  results  for 
lines  B-F,  indicating  no  significant  subsurface  anomalous  condi¬ 
tions.  A  synopsis  of  features  of  the  GPR  records  for  lines  C-F 
is  given  in  Table  2.  Three  interpretation  guidelines  used  in 
preparing  this  synopsis  are: 

a.  Interfaces  between  different  materials  will  produce 
reflection  events  on  the  graphic  record  which  vary 
in  record  time  in  relation  to  variations  in  depth 
to  the  interface  (  depth  -  0.5  x  t  x  C  ;  where  t  is 
the  event  time  from  the  record,  and  C  is  the  speed 
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Figure  n.  QPR  record  for  line  E. 


Table  2 


Synopsis  of  GBR  Results  Lines  C-F 


LXMB  C 

Anomaly  Looatlon 


Description 


General  Thin  surface  layer  (1  to  3  ft 

thick)  along  entire  line. 

Below  this  reflector,  mater¬ 
ials  are  generally  uniform, 
with  occasional  diffractions 
and  small  areal  extent  reflec¬ 
tion  events. 


C2-C6 


C49-C63 

C1-,C11,C36, 

C42,C60 


Disturbed  area/shallow  depres¬ 
sion,  approx.  5  ft  depth. 

Shallow  subsurface  depression. 

Small,  shallow  dlffractors, 
3-10  depth 


LIME  D 

Anomaly  Location 


Description 


General 


D31-D35 

D9 


Thin  surface  layer  (1  to  3  ft 
thick)  along  entire  line. 

Very  uniform  conditions  to 
depth  of  approx.  20  ft. 

shallow  disturbed  area, 

3-5  ft  depth. 

Small,  shallow  diffractors. 


Anomaly  Location 
General 

E1-E13 

E16,E30 


LIME  E 

Description 

Thin  surface  layer  (1  to  3  ft 
thick)  along  line. 

Presence  of  pavement  produces 
ringing  on  record. 

Small,  shallow  diffractors. 

(Continued) 
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Table 

Anomaly  Location 
General 

F2 

F13-F21 

F14 

F42 

F65 


2  (Concluded) 

LIMB  F 

Description 

Thin  surface  layer  (1  to  4  ft 
thick)  along  entire  line  and 
generally  uniform  conditions 
to  depth  of  approx.  20  ft. 

Several,  small,  closely-* 
spaced  dlff rectors,  approx. 

5  ft  depth. 

Ringing  on  record  due  to  sur¬ 
face  pavement. 

Small  dlffractors,  8-10  ft 
depth . 

Small  diffraction  apparently 
due  to  water  well  casing, 
located  approx.  15  ft  horlz. 

Small  depression,  approx. 

5  ft  depth,  with  apparent 
diffraction  below. 
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of  alectronagnetic  waves  in  the  subsurface  or 
approximately  4  x  10®  ft/s  or  1.2  x  lo®  m/s); 

b.  Changes  in  material  type  or  condition  (e.g.,  water 
content,  grain  size,  cementation,  clay  content) 
along  the  interface  will  produce  changes  in  the 
character/appearance  of  the  reflection  event  and 

in  some  cases  a  breakup  in  the  continuity  of  the 
event ; 

c.  Localized  features,  e.g.  cavities,  boulders/cobbles, 
caliche  zones,  low  density  zones,  high  porosity/ 
water  content  zones,  etc.,  with  size  of  the  order 

of  the  electromagnetic  wavelength  (approx.  3  ft) 
will  produce  a  hyperbolic-shaped  event  (diffraction) 
on  the  record  which  opens  downward  in  time. 

29.  If  the  interpretation  of  the  GPR  records  for  lines  C-P 
is  simple  due  to  lack  of  detail  on  the  records,  then  the  inter¬ 
pretation  of  the  records  for  the  sinkhole  GPR  surveys  is  very 
difficult  due  to  the  considerable  detail  exhibited  in  the 
records.  Figure  12  is  the  GPR  record  for  E-W  line  10,  located 
50  ft  N  of  the  sinkhole.  This  record  illustrates  the  detail  and 
complexity  exhibited  by  most  of  the  sinkhole  GPR  records.  The 
most  prominent  feature  of  Figure  12  is  a  basin-like  structure 
centered  on  N-S  cross  line  4;  depth  to  the  apparent  bottom  of  the 
structure  is  6  ft.  Considering  only  the  structure  revealed  in 
Figure  12,  it  could  be  a  filled  sinkhole  or  shallow  channel. 
Beneath  the  basin  structure  is  a  complex  of  diffractions  and  dip¬ 
ping  events,  and  the  trough-like  depression  between  cross  lines  3 
and  4  is  possible  evidence  of  a  migration/piping  pathway.  Two 
other  features  between  cross  lines  4  and  5  resemble  vertical 
pipes.  The  GPR  record  for  N-S  line  5  is  presented  in  Figure  13, 
and  exhibits  a  basin-like  structure  which  extends  S  of  the  extent 
of  the  GPR  survey  line.  Also,  the  vertical  pipe  feature  nearest 
cross  line  5  in  Figure  12  appears  in  Figure  13  near  cross  line  10 
as  expected  and  suggests  that  the  feature  is  'quasi-circular'  in 
horizontal  section  and  not  elongate. 

30.  In  order  to  summarize  the  sinkhole  GPR  surveys,  a 
three  dimensional  fence  diagram  of  the  major  structures  revealed 
by  the  records  is  shown  in  Figure  14,  and  a  plan  map  of  smaller 
scale  anomalies,  such  as  the  vertical  pipe  features  discussed 
above,  is  presented  in  Figure  15.  The  trends  of  the  major  struc¬ 
ture  in  Figure  14  are  (l)a  dip  to  the  east  produced  by  2-3  N-S 
'escarpment- like'  steps  and  (2) a  more  uniform  dip  to  the  south. 
The  overall  dip  of  the  major  reflector  is  to  the  SE.  The  survey 
area  is  too  small  to  determine  if  the  major  structure  is  a  buried 
channel  trending  to  the  SE  or  a  large,  filled  sinkhole.  Based  on 
elevations  of  the  top  of  the  conglomerate  (1180  ft  in  RDH-1)  in 
drillholes  located  several  hundred  feet  to  the  W  of  the  sinkhole, 
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Figure 


CROSS  LINES 


Figure  13.  GPR  sinkhole  survey  line 
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Figure  15.  Map  of  localized  anomalous  features  from  sinkhole  GPR  surveys. 
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the  reflecting  basement  In  the  sinkhole  6PR  surveys  which  lies 
5-10  ft  below  the  surface  (nominal  surface  elevation  —  1190  ft) 
is  likely  the  conglomerate.  Two  types  of  anomalous  features  are 
located  in  Figure  15,  hyperbolic  radar  signatures  (possible 
cavities)  and  signatures  of  possible  vertical  pipes/ col lapse  fea¬ 
tures  discussed  above  for  lines  5  a)  d  10.  There  are  two  align¬ 
ments  of  anomalies  which  may  be  significant:  (1)  an  alignment  of 
features  between  the  N-S  survey  lines  3  and  4,  which  includes  the 
existing  sinkhole;  (2)  a  possible  trend  of  anomalies  from  NE  to 
SW,  indicated  by  the  dashed  line  in  Figure  15. 

Electromagnetic  Conductivity 

31.  EM  surveys  were  conducted  along  lines  C,E,  and  F. 
Results  of  the  EM  surveys  are  given  in  Figures  16-21;  the  first 
of  the  two  figures  for  each  line  includes  the  three  HD  profiles, 
while  the  second  includes  the  three  VO  profiles.  For  line  C, 
high  conductivity  anomalies  exist  at  each  end  of  the  line,  with 
nearly  constant/uniform  conditions  vertically  and  horizontally  in 
the  central  portion  of  the  line.  The  anomalies  at  each  end  of 
line  C  extend  beyond  the  limits  of  the  survey.  Using  the  ap¬ 
proximate  depths  of  investigation  of  the  EM  method  discussed  in 
Paragraph  19,  a  cross  section  of  the  EM  anomalies  is  shown  in 
Figure  22.  The  dashed  lines  which  are  parallel  to  the  surface 
indicate  the  approximate  depth  of  investigation  for  the  three 
vertical  dipole  coil  spacings.  Also  shown  in  the  cross  section 
are  the  locations  of  the  exploratory  borings,  which  were  con¬ 
ducted  during  and  shortly  after  the  geophysical  field  program, 
and  the  ends  of  the  concrete  cutoff  wall  and  grout  curtain. 
Within  the  large  anomalous  areas  Cl  and  C2,  conductivities  are  in 
excess  of  100  mmhos/m  (10  ohm-m  resistivity) .  Outside  the 
anomalous  areas,  the  conductivity  is  15-20  mmhos/m  (50-65  ohm-m 
resistivity),  with  minor  anomalies  as  high  as  30  mmhos/m  (33 
ohm-m  resistivity) .  The  Cl  and  C2  anomalies  are  most  likely  due 
to  high  water  content,  although  high  clay  content  is  also  a  pos¬ 
sibility. 


32.  The  exploratory  borings  encounter  gravel  and  cobbles, 
which  is  interpreted  to  be  the  conglomerate,  at  a  nearly  constant 
elevation  along  the  survey  line.  Above  the  conglomerate,  the 
borings  encounter  silt,  which  varies  from  dry  to  moist,  and 
caliche,  which  varies  from  scattered  small  nodules  to  zones 
several  inches  thick.  Borings  U-5  and  U-10  were  anomalous  in 
that  U-5  encountered  wet  conditions  at  the  top  of  the  con¬ 
glomerate  and  U-10  encountered  6  ft  of  wet  silt  above  the  con¬ 
glomerate.  Borings  U-9  and  U-10  were  placed  to  Investigate 
anomaly  C2.  The  wet  zone  in  U-10  supports  the  interpretation  of 
a  high  water  content  explanation  for  anomalies  Cl  and  C2. 
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Figure  16.  EM  conductivity  profiles,  line  C,  horizonta 
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conductivity  profiles,  line  C,  vertical  dipoles. 


Figure  1?.  EM  conductivity  profiles,  line  E,  horizontal  dipoles. 
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Figure  21.  EM  conductivity  profiles,  line  F,  vertical  dipoles. 


CONDUCTIVITY  ANOMALIES 


Figure  22.  EM  anomaly  cross  section,  line 


33.  For  line  E,  a  high  conductivity  anomaly  is  indicated 
from  the  beginning  of  the  line  to  approximately  the  220  ft  posi¬ 
tion.  Beyond  this  anomaly,  the  conductivity,  except  for  the 
40  m  VD  data,  is  approximately  the  same  as  in  the  central  region 
of  line  C.  The  40  m  vertical  dipole  data  in  the  220-350  ft 
region  is  elevated  by  nearly  50  percent  relative  to  the  rest  of 
the  EM  data,  which  could  Indicate  higher  water  contents  at  depths 
greater  than  about  100  ft. 

34.  The  EM  data  for  line  F,  Figures  20  and  21  is  highly 
disturbed  by  cultural  features:  existing  toilet  facilities  at 
location  450-475  ft;  a  deep  well  with  metal  casing  at  location 
410  ft;  an  unmarked,  buried  cultural  feature  at  location  130  ft 
(possibly  a  pipe  associated  with  previous  toilet  facilities  near 
this  location) .  Along  the  segments  of  the  survey  line  away  from 
the  cultural  feature  influence,  the  conductivity  is  uniform  and 
approximately  equal  to  the  background  levels  for  lines  C  and  E. 
A  shallow  anomalous  feature  (<20  ft  depth)  exists  at  the  right 
end  of  the  survey  line  (location  600-700  ft),  and  a  deeper 
anomalous  feature  is  indicated  at  the  left  end  of  the  line 
(location  0-100  ft)  .  Plan  locations  of  all  the  major  EM 
anomalies  are  indicated  on  Figure  23. 

Electrical  Resistivity 

33.  Data  from  the  individual  soundings  for  the  pole-dipole 
resistivity  survey  along  line  c  are  presented  in  Appendix  C.  The 
first  set  of  sounding  plots  in  the  appendix  are  'raw'  apparent 
resistivity  values  plotted  at  the  midpoint  of  the  potential 
electrode  pair  locations,  while  the  second  set  of  data  results 
from  a  center-weighted,  three  point  moving  filter  applied  to  the 
raw  data  to  produce  smoothed  sounding  curves.  Anomalies  are 
identified  relative  to  the  smoothed  curves.  Using  the  graphical 
interpretation  procedure  discussed  in  Paragraph  21,  the  resis¬ 
tivity  anomaly  cross  section  shown  in  Figure  24  results.  It  is 
highly  significant  that  anomalies  HI  and  R2  correlate  both  in 
sense  (high/low)  and  location  to  the  EM  conductivity  anomalies  Cl 
and  C2.  The  small,  low  resistivity  anomaly  R3  is  apparently  con¬ 
firmed  by  boring  U-5,  which  encountered  wet  conditions  at  the  top 
of  the  conglomerate.  Of  anomalies  R3-R7,  only  R5  is  possibly  ex¬ 
pressed  in  the  EM  conductivity  data.  The  high  resistivity 
anomalies  are  likely  caused  by  dryer  conditions  in  the  silt 
and/or  zones  Wfith  higher  caliche  concentration.  The  EM  method 
does  not  have  the  vertical  resolution  of  the  pole-dipole  method 
and  it  is  likely  that  the  high  resistivity  anomalies  are  masked 
in  the  EM  data  by  the  large  response  to  anomc  ;.es  Cl  and  C2. 
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Figure  23.  EM  anomaly  locations. 
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RESISTIVITY  ANOMALIES 


Figure  24.  Resistivity  anomaly  cross  section,  line  C. 


34.  The  Bouguer  gravity  results  for  lines  are  shown  in 
Figures  25-29  respectively.  The  results  are  relative  to  the  base 
station  for  line  C  (C37  or  360  ft);  i.e.,  the  data  for  each  line 
are  corrected  to  the  latitude  and  elevation  for  C37,  the 
reference  gravity  value  for  C37  is  subtracted  from  all  the  data, 
and  then  the  data  are  terrain  corrected.  Beference  or  'best 
fit'  curves  through  the  data  are  given  for  each  line.  The  data 
are  scattered  relative  to  these  curves,  representing  actual 
anomalies  caused  by  small,  shallow  subsurface  features  (typically 
less  than  20  ft  or  6  m  deep),  data  measurement  errors,  and  in 
many  cases  errors  due  to  incomplete  terrain  correction.  Incom¬ 
plete  terrain  corrections  in  this  case  refers  to  incomplete 
definition  of  terrain  variations  close  to  each  measurement  loca¬ 
tion,  say  within  a  10  ft  (3m)  radius.  Incomplete  terrain  cor¬ 
rection  errors  will  affect  the  data  for  lines  C  and  D  more  than 
for  lines  D,  E,  and  F.  The  complete  set  of  measured  and 
processed  gravity  data  is  included  in  Appendix  D  in  the  form  of 
computer  printouts. 

35.  The  anomalies  caused  by  the  curves  in  Figures  25-29 
are  caused  by  large  scale  and/or  deep-seated  features  (relative 
to  the  small  anomalies  discussed  above) .  The  old  addage,  that 
anomalies  inevitably  occur  near  the  ends  of  survey  lines,  is 
dramatically  illustrated  by  the  gravity  data  for  lines  B,  c,  and 
E,  where  low  gravity  anomalies  occur  near  the  ends  of  the  lines. 
Low  gravity  anomalies  also  appear  in  the  data  for  lines  D  and  F, 
but  these  lows  are  more  completely  defined.  Figure  30  is  a  plan 
map  showing  the  spatial  extent  of  the  low  gravity  anomalies.  The 
center  of  the  anomaly  can  be  defined  for  lines  D  and  F  and  for 
one  of  the  anomalies  on  line  C.  These  anomaly  centers  are  indi¬ 
cated  on  Figure  30;  and  the  double  dashed  line  indicates  a  sig¬ 
nificant  alignment  of  low  gravity  anomalies  (Trend  II) ,  possibly 
caused  by  a  low  density  anomalous  feature  extending  through  the 
right  abutment  and  passing  underneath  the  dam.  Another  possible 
alignment  of  low  gravity  anomalies  is  indicated  in  Figure  30  near 
the  right  ends  of  lines  B,  C,  and  E  (Trend  1),  also  apparently 
caused  by  a  feature  extending  underneath  the  dam. 

36.  The  low  gravity  anomalies,  for  which  centers  were  in¬ 
dicated  above,  are  complete  enough  for  depth  estimates  to  be 
made.  Depth  estimates  are  made  for  a  particular  assumption  of 
feature  geometry.  The  anomalous  feature  is  assumed  to  be  elon¬ 
gated,  with  a  compact,  constant  cross  sectional  area  and  to  have 
a  sub-horizontal  orientation,  e.g.,  a  long,  nearly  horizontal 
circular  cylinder.  Using  this  assumption,  the  following  depth 
estimates  (to  center  of  the  anomalous  feature)  result: 
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Terrain  corrected  Bouguer  gravity  profiie.  line  C 


LINE  D 
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Figure  27.  Terrain  corrected  Bouguer  gravity  profile,  line 
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Figure  28.  Terrain  corrected  Bouguer  gravity  prcfiie,  line 


LINE  F 


errain  corrected  Bouguer  gravity  profile,  line 


Line 


Anomaly 

Center 


Anomaly  Magnitude 
(microgals) 


Depth 

(ft) 


c 

C71 

(700 

ft) 

-  110 

190 

F 

F43 

(420 

ft) 

-  105 

140 

D 

D47 

(460 

ft) 

-  65 

100 

The  depth  estimates  are  arranged  according  to  decreasir<g  depth. 
Since  the  anomaly  magnitude  apparently  increases  with  increasing 
depth,  the  anomalous  feature  must  be  larger  in  physical  size 
and/or  density  contrast  as  it  passes  successively  under  lines  D, 
F,  C  and  B. 

37.  The  above  interpretation  based  on  the  simple  assump-* 
tion  of  a  model  geometry  is  not  unique.  First,  for  the  computed 
depth,  in  each  case  only  the  product  of  density  contrast  and  mean 
radius  can  be  computed  (Butler  1980) ,  where  density  contrast 
refers  to  the  difference  between  the  density  of  the  anomalous 
feature  and  the  density  of  the  surrounding  material.  The  follow¬ 
ing  tabulations  are  calculations  of  the  mean  radius  of  the 
anomalous  feature  for  selected  density  contrasts  for  two  of  the 
cases  presented  above: 

Line  C  —  Depth  to  Center  190  ft  t58  m) 

Density  Contrast  Mean  Radius 
g/cm^  ft  (rc) 


-  0.1 

130 

(40) 

-  0.2 

90 

(27) 

-  0.3 

74 

(23) 

-  0.4 

64 

(20) 

-  0.5 

57 

(17) 

-  2.2 

27 

(8) 

Line  D  —  Depth  to  Center  =  100  ft  (30  m) 

Density  Contrast  Mean  Radius 

g/cTxr  ft  (m) 


71  (22) 
50  (15) 
41  (13) 
36  (11) 
32  (10) 
15  (5) 


-  0.1 
-  0.2 

-  0.3 

-  0.4 

-  0.5 

-  2.2 
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CONSTRUCTION  9ASELINE 


The  -*  2*2  g/ota?  density  contrast  is  Included  as  an  extreme  case, 
representing  an  air-filled  cavity  in  the  conglomerate;  a  water- 
filled  cavity  would  have  a  density  contrast  of  -  1.2  q/cvr ,  A 
more  likely  density  contrast  of  uncemented  relative  to  cemented 
conglomerate  or  a  zone  in  which  piping  has  occurred  is  -0.2  to 
-0.3  g/cvr,  A  second  way  in  which  the  interpretation  may  not  be 
unique  is  that  the  assumption  of  feature  geometry  may  be  invalid, 
other  assumptions  of  geometry  will  result  in  shallower  depths 
(Butler  1980) . 

38.  A  Bouguer  anomaly  contour  map  for  the  sinkhole 
microgravity  survey  is  shown  in  Figure  31,  where  the  circled  num¬ 
bers  correspond  to  the  GPR  survey  lines.  An  anomaly  caused  by 
the  small  existing  sinkhole  does  not  appear  on  the  gravity  map 
due  to  insufficient  gravity  measurements  in  the  immediate 
vicinity  of  the  sinkhole.  A  three  dimensional  representation  of 
the  Bouguer  anomaly  data  is  presented  in  Figure  32  from  the  same 
perspective  as  the  GPR  fence  diagram  in  Figure  14.  The  gravity 
data  surface  has  the  same  trend  as  the  major  GPR  reflector,  an 
overall  dip  to  the  SE  (low  gravity  anomaly  to  the  SE) .  The  low 
gravity  trend  to  the  SE  supports  the  suggestion  that  the  major 
GPR  reflector  dipping  to  the  SE  is  the  conglomerate,  since  the 
conglomerate  should  represent  a  positive  density  contrast  rela¬ 
tive  to  the  silt.  A  localized  low  gravity  anoraaly  near  location 
(30,30)  is  well  defined  and  correlates  well  in  location  with  a 
hyperbolic  GPR  feature  (cavity;  and  the  NE-SW  significant  GPR 
anomaly  trend  noted  earlier. 


Integrated  Assessment 

39.  An  integrated,  complementary  geophysical  program  was 
conducted  at  the  site,  however  the  results  of  the  individual 
methods  varied  greatly  in  their  applicability  and  contribution  to 
the  program.  The  shallow,  seismic  reflection  method  was  least 
successful  at  the  site  and  contributed  least  to  the  program. 
Lack  of  identifiable  silt/congloroerate  reflections  in  the  results 
of  the  seismic  tests  likely  indicates  that  the  conglomerate  is 
uncemented  in  the  three  locations,  which  is  consistent  with  the 
gravity  survey  results  which  have  low  gravity  (low  density) 
anomalies  in  the  three  locations. 

40.  The, GPR  has  a  maximum  effective  depth  of  investigation 
at  the  site  of  approximately  30  ft.  GPR  records  along  the  toe  of 
the  dam  and  on  the  right  abutment  are  relatively  featureless, 
i.e.,  the  records  indicate  relative  uniform  conditions  laterally 
and  vertically  to  a  depth  of  30  ft.  In  contrast,  the  GPR  results 
In  the  vicinity  of  the  sinkhole  in  the  floor  of  the  reservoir 
were  remarkable  in  the  amount  of  detail  exhibited  in  the  records. 
The  sinkhole  GPR  records  indicate  a  prominent  reflector  dipping 
CO  the  SE  and  isolated  vertical  pipe  features  and  possible 


53 


DfSTANCE.  FT 


/  oBojo]  w'Aj}  AOjg 


Figure  32.  Three-dimensional  representation  of  gravity  data  around  sinkhole. 
Circled  numbers  refer  to  the  GPR  survey  lines. 
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oavltlss  (Figures  14  and  IS) .  The  sllt/conglomerate  interface  is 
a  viable  candidate  for  the  reflector  in  the  vicinity  of  the 
sinkhole.  7he  results  and  interpretation  of  the  sinkhole  @PR 
surveys  and  consistent  with  the  sinkhole  microgravlty  survey 
results  and  Interpretation  (Figure  32). 

41.  Along  the  upstream  toe  of  the  dam  in  the  right  abut¬ 
ment  area,  line  C,  results  and  Interpretations  of  the  electromag¬ 
netic  conductivity,  electrical  resistivity,  and  microgravity 
surveys  are  quite  consistent.  All  three  methods  indicate 
anomalous  conditions  at  each  end  of  line  C.  The  subsurface  fea¬ 
tures  causing  the  anomalies  are  interpreted  as  low  density,  high 
water  content  zones.  The  EM  and  ER  interpretations  indicate  that 
the  anomalous  features  extend  from  depths  of  50-100  ft  to  depths 
of  150-200  ft  or  greater  (Figures  22  and  24) .  one  interpretation 
of  the  microgravity  results,  which  assumes  that  the  anomaly  at 
the  north  end  of  line  C  is  caused  by  a  cylindrical  cross-section 
feature,  gives  a  depth  to  anomaly  center  of  190  ft  and  a  mean 
anomaly  radius  of  74  ft  (90  ft)  for  an  assumed  density  contrast 
of  -  0.3  g/cxer  (-  0.2  g/cm^).  This  microgravity  interpretation 
is  consistent  with  the  EM  and  ER  interpretations.  The  ER  and  EM 
interpretations  are  confirmed  in  two  locations  along  line  C  by 
results  of  the  exploratory  drilling  program. 

42.  The  microgravity  survey  results  and  Interpretations 
indicate  two  anomalous  trends  which  apparently  pass  under  the  dam 
(Figure  30,  Trends  I  and  II).  The  south  cnmost  trend.  Trend  1, 
is  consistent  with  the  possible  trend  in  bleated  by  EM  anomalies 
shown  in  Figure  23. 

43.  The  locations  of  the  anomalous  zones  in  relation  to 
the  right  abutment  terminations  of  the  concrete  cuteff  wal'.  and 
grout  curtain  is  likely  not  coincidental,  and  strongly  suc.gests 
that  the  cause  of  the  anomalies  is  seepage  related.  A  seepage- 
related  cause  for  the  anomalous  zones  is  supported  by  the  results 
of  the  self  potential  (SP)  surveys  and  other  geotechnical  inves¬ 
tigations  conducted  during  the  1984  test  fill  of  the  reservoir 
(Butler,  Wahl  and  Sharp  1984).  Two  long  self  potential  survey 
lines  were  established  which  included  the  present  lines  C  and  F. 
For  the  SP  line  which  paralled  the  dam  centerline  along  the 
cutoff  wall/grout  curtain  alignment  (including  the  present  line 
C) ,  anomalous  seepage  was  indicated  near  the  dam/right  abutment 
contact  and  termination  of  the  cutoff  wall.  Along  the  SP  line 
perpendicular  to  the  dam  (including  the  present  line  F)  , 
anomalous  seepage  was  indicated  very  close  to  the  center  of  the 
low  gravity  finomaly  on  line  F. 
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PART  IV  :  CONCLUSIONS  AND  RECOMMENDATIONS 


ConcluslonG; 

43.  An  integrated,  complementary  geophysical  program 
directed  to  the  detection  and  delineation  of  anomalous  conditions 
in  the  right  abutment  of  Mill  Creek  Dam,  Walla  Walla,  Washington, 
was  planned  and  conducted.  During  and  shortly  after  the 
geophysical  field  program,  exploratory  drilling  was  conducted 
along  the  upstream  toe  of  the  dam  in  the  right  abutment  area. 
The  following  major  conclusions  result  from  the  integrated  as¬ 
sessment  of  the  geophysical  program  and  exploratory  drilling: 

a.  Two  anomalous  trends  apparently  extend  from  the 
reservoir  through  the  right  abutment  and  under  the 
dam;  interpreted  depths  of  the  two  anomalous  feature 
trends  place  the  features  likely  within  the  con¬ 
glomerate;  the  interpreted  features  have  low  density 
and  high  water  content;  one  of  the  anomalous  trend 
interpretations  is  confirmed  by  an  exploratory 
borehole  which  encountered  a  5  ft  wet  zone  in  the 
silt  just  above  the  conglomerate; 

b.  Surveys  around  a  small  sinkhole  in  the  reservoir 
floor  reveal  a  conglomerate  surface  which  dips  to 
the  S-SE;  several  vertical  pipe  features  and  small 
cavities  are  detected  within  the  conglomerate, 
consistent  with  a  scenario  of  vertical  seepage  and 
piping  of  silt  into  the  conglomerate; 

3.  The  shallow  subsurface  in  the  right  abutment  area  to 
a  depth  of  30  ft  is  relatively  uniform  laterally; 
the  exploratory  drilling  confirms  the  geophysical 
interpretation . 

44.  Examination  of  all  the  results,  particularly  along 
line  C,  suggests  that  the  anomalous  conditions  and  trends  are 
controlled  or  at  least  strongly  influenced  by  the  cutoff  wall  and 
grout  curtain.  A  seepage-related  cause  for  the  anomalous  trends 
is  strongly  suggested  by  the  present  results  and  supported  by 
results  of  self  potential  surveys  and  other  geotechnical  inves¬ 
tigations  conducted  during  the  1984  reservoir  test  fill. 


Recommendations 

45.  The  conclusions  (1)  that  two  low  density,  high 
water  content  anomalous  zones  pass  through  the  right  abutment  and 
underneath  the  dam  and  (2)  that  these  anomalous  zones  are  caused 
by  seepage  related  processes  are  interpretations,  based  on  exten- 
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slvo  geophysical  and  geotechnical  data^  and  are  considered  well 
founded.  However,  the  interpretations  were  only  marginally  con¬ 
firmed  by  direct  sampling  of  material  during  the  exploratory 
drilling  program.  If  further  exploratory/confirmatory  drilling 
Is  undertaken,  a  minimum  recommended  program  Is  presented  in 
figure  33,  where  the  basis  for  the  recommendation  Is  Indicated. 
The  drilling  program  Ideally  should  Include  sampling  Immediately 
above  and  within  the  conglomerate  and  holes  should  extend  to  the 
top  of  the  basalt.  Drilling  In  the  conglomerate,  however,  will 
be  difficult  and  expensive  and  the  results  possibly  ambiguous  due 
to  the  drilling  and  sampling  difficulties.  Selective,  shallow 
drilling  In  the  vicinity  of  the  sinkhole  to  depths  of  at  least  30 
ft  Is  recommended  would  be  more  easily  achievable. 

46.  Depending  on  project  fate  decisions  to  be  made  by  the 
District  and  Division,  future  geophysical  investigations  may  be 
justified.  Using  the  present  geophysical  data  as  a  baseline,  fu¬ 
ture  surveys  could  detect  changes  In  subsurface  conditions  in  the 
right  abutment  area. 
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Figure  33.  Recommpnded  minimum  confirmatory  drilling  program 


REFERENCES 


Ballard,  Robert  F.  1983.  “Cavity  Detection  and  Delineation 
Research,  Report  5,  Electromagnetic  (Radar)  Techniques  Applied  to 
Cavity  Detection,"  Technical  Report  GL-83-1,  US  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS. 

Bates,  E.R.  1973.  "Detection  of  Subsurface  Cavities,"  Miscel¬ 
laneous  Paper  S-73-40,  US  Army  Engineer  Waterways  Experiment  Sta¬ 
tion,  Vicksburg,  MS. 

Bristow,  C.M.  1966.  "A  New  Graphical  Resistivity  Technique  for 
Detection  of  Air-Filled  Cavities,"  Studies  in  Speleology,  v.  7, 
pp  204-227. 

Butler,  Dwain  K.  1980.  "Microgravimetric  Techniques  for 
Geotechnical  Applications,"  Miscellaneous  Paper  GL-80-13,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Butler,  Dwain  K.  1986.  "Military  Hydrology,  Report  10,  Assess¬ 
ment  and  Field  Examples  of  Continuous  Wave  Electromagnetic  Sur¬ 
veying  for  Ground  Water,"  Miscellaneous  Paper  EL-79-6,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Butler,  Dwain  K. ,  and  Yule,  Donald  E.  1984.  "Microgravity  Sur¬ 
vey  of  Wilson  Dam  Powerplant  Switchyards,  Florence,  Alabama," 
Miscellaneous  Paper  GL-84-2,  US  Army  Engineer  Waterways  Experi¬ 
ment  Station,  Vicksburg,  MS. 

Butler,  D.K.,  Gangi,  A.F.,  Wahl,  R.E.,  Yule,  D.E.,  and 
Barnes,  D.E.  1982.  "Analytical  and  Data  Processing  Techniques 
for  Interpretation  of  Geophysical  Survey  Data  with  Special  Ap¬ 
plication  to  Cavity  Detection,"  Miscellaneous  Paper  GL-82-16,  US 
Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Butler,  Dwain  K. ,  Wahl,  Ronald  E. ,  and  Sharp,  Michael  K.  1984. 
"Geophysical  Seepage  Detection  Studies,  Mill  Creek  Dam,  Walla 
Walla,  Washington,"  Miscellaneous  Paper  GL-84-16,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Department  of  the  Army  1979.  "Geophysical  Exploration,"  Engineer 
Manual  EM  1110-1-1802,  Office  of  the  Chief  of  Engineers, 
Washington,  D.C. 

Dobecki,  T.L.,  Mueller,  T.L.,  Savage,  M.B.  1989.  "High- 
Resolution  Seismic  Reflection  Investigations,  Beaver  Dam, 
Arkansas,"  Technical  Report  REMR-GT-10,  US  Army  Engineer  Water¬ 
ways  Experiment  Station,  Vicksburg,  MS. 

Telford,  W.M. ,  Geldart,  L.P.,  Sheriff,  R.E.,  and  Keys,  D.A. 
1976.  Applied  Geophysics.  Cambridge  University  Press,  New  York. 


60 


Walla  Walla  District,  Corps  of  Engineers.  1987.  "Mill  Creek 
Lake,  Washington,  General  Design  Memorandum  No.  5,  Project 
Rehabilitation:  Supplement  No.  2,  Storage  Dam  Seepage  Control," 
U.S.  Army  Engineer  District,  Walla  Walla,  Walla  Walla,  WA. 


APPENDIX  A 

SEISMIC  REFLECTION  SURVEY,  MILL  CREEK  DAM 

Great  Plains  Geophysical 
Lawrence,  Kansas 


SEISMIC  REFLECTION  SURVEY 
FEASIBILITY  STUDY 


AT 

MILL  CREEK  DAM 
NEAR 

WALLA  WALLA,  WASHINGTON 


A  REPORT  PREPARED  FOR 
THE  US  ARMY  CORPS  OF  ENGINEERS 
WATERWAYS  EXPERIMENT  STATION 
VICKSBURG,  MISSISSIPPI 


BY 

GREAT  PLAINS  GEOPHYSICAL,  INC. 
2913  WESTDALE  ROAD 
LAWRENCE.  KANSAS  66049 


JULY  10,  1989 


A3 


Seismic  reflection  was  used  in  conjunction  with  a  variety  of 
other  geophysical  techniques  at  Mill  Creek  Dam  near  Walla  Walla, 
Washington  to  assess  the  embankment  and  foundation  conditions  in 
the  dam/right  abutment  contact  area.  The  reflection  survey  was 
designed  to  detect  the  suspected  cavities  in  a  silt  deposit  (the  silt 
deposit  is  the  foundation  of  the  dam  as  well  as  the  right  abutment. 
The  shallow  stratigraphy  of  the  study  area  includes  a  surface  silt 
layer  varying  in  thickness  from  40  to  100  ft  overlying  an 
approximately  100  ft  thick  conglomerate  which  has  a  basalt  at  its 
basal  contact.  The  integrity  of  the  dam  could  be  threatened  if 
sinkholes  similar  to  those  detected  in  the  floor  of  the  reservoir 
continue  to  develop.  The  mechanism  responsible  for  the  cavity 
formation  is  uncertain.  It  is  suspected  that  piping  of  silt  into  more 
pervious  zones  of  the  conglomerate  could  be  responsible.  The 
conglomerate  may  also  possess  cavities  as  a  result  of  piping  of  silts 
from  within  pockets  in  the  conglomerate. 

The  data  were  recorded  on  a  Bison  9024  digital  floating  point 
seismograph.  The  seismograph  has  the  capability  to  analog  filter 
with  user  selectable  high  and  low  cut  filters  on  4  hz  increments.  The 
seismograph  generates  a  16  bit  digital  word  with  a  two  step  floating 
point  amplifier  (0  or  30  dB)  and  a  16  bit  A/D  converter. 

Three  different  sources  were  tested  at  several  different  filter 
settings  at  three  locations  (shots  1-57  at  site  #1,  shots  58-81  at  site 
#2,  shots  83-101  at  site  #3)  on  the  north  end  of  the  reservoir,  for  a 
total  of  101  test  shots.  There  are  no  obvious  reflectors  at  any  of  the 
three  sites  with  any  of  the  parameters  used.  No  reflection  was 
observed  from  the  top  of  the  conglomerate  or  from  within  the 
conglomerate. 

The  20  pound  sledge  hammer,  the  downhole  30.06  rifle,  and 
the  downhole  .50  caliber  rifle  were  all  individually  used  as  energy 
sources.  The  body  wave  energy  from  the  rifle  sources  has  a  slightly 
higher  frequency  with  less  ground  roll,  although  the  hammer  data 
are  a  little  less  'ringy'.  There  was  little  or  no  seismic  energy  above 
300  Hz,  so  the  use  of  40  Hz  versus  100  Hz  geophones  was  not  critical, 
except  to  the  extent  that  the  3  L-28E  40  Hz  geophones  (which  were 
the  only  phones  tested)  have  higher  voltage  output  than  the  single  L- 
40A  100  Hz  geophones,  given  the  same  input  ground  velocity.  Lo- 
cut  filters  in  the  200  Hz  range  were  efficient  in  removing  a 
significant  amount  of  ground  roll. 
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Analysis  of  first  arrival  information  (direct  wave  and  refraction 
velocities  versus  offset)  suggest  the  presence  of  a  shallow  refractor 
at  a  depth  of  10  feet  or  less.  The  P-wave  velocity  in  the  surface 
layer  (soil)  is  in  the  600  to  1000  foot/second  range.  The  velocity  in 
the  refracting  layer  (loess  and/or  silt)  is  in  the  1400  to  3200 
foot/second  range.  Maximum  offsets  were  not  sufficient  to 
determine  the  refractive  thickness  of  the  loess  layer  or  the  velocity 
of  the  underlying  conglomerate  or  basalt. 

Two  events  with  some  subtle  reflection  characteristics  are 
identifiable  at  two  different  sites  and  at  two  different  times.  The 
shallower  of  the  two  possible  reflections  occurs  at  site#l  at  a  time  of 
approximately  100  msec.  This  event  is  suspect  of  not  being  a  true 
reflection  due  to  the  highly  irregular  first  arrival  on  walkaways  at 
this  site  (probably  the  result  of  the  slurry  wall)  and  the  lack  of  good 
coherency  for  more  than  about  10  traces.  This  shallower  event 
shows  a  zero-distance  time-intercept  of  92  msec,  a  depth  of  171  feet, 
and  a  normal-moveout  velocity  of  3730  ft/sec.  The  correlation 
coefficient  of  a  least-squares  hyperbola  fit  to  the  data  is  0.998. 
Generally,  anything  above  0.99  is  considered  to  be  a  hopeful  sign  of  a 
true  reflector.  The  first  57  test  shots  were  fired  at  test  site  #1  which 
has  the  slurry  wall  present  as  a  complication.  As  a  result,  skepticism 
exists  concerning  results  at  site  #1. 

The  suggestion  of  reflection  information  can  be  interpreted  on 
walkaways  from  site  #3  (with  a  200  hz  or  higher  lowcut  filter)  at  a 
time  of  approximately  190  msec.  The  event  has  coherency  across 
traces  with  offsets  less  than  about  80  ft.  The  small  source  to  receiver 
offset  of  the  two  dozen  or  so  traces  with  coherency  results  in  no 
detectable  moveout.  This  potential  reflector  has  a  very  "broken  up" 
look  which  if  real  must  be  an  indication  of  significant  near-surface 
static  anomalies.  The  nearly  flat  appearance  of  the  event  restricts 
the  effectiveness  of  an  NMO  calculation  to  determine  depth,  origin 
time,  velocity,  and  correlation  to  a  hyperbola.  The  frequency  content 
and  sporadic  nature  of  the  coherency  suggest  possible  spatial 
aliasing.  The  only  true  way  to  determine  if  this  event  is  the  result  of 
spatial  aliasing  is  by  moving  the  source  to  a  slightly  different  offset 
during  acquisition.  This,  of  course  is  not  possible  at  this  time. 

We  see  no  sign  of  reflections  at  site  #2.  It  should  be  pointed 
out,  however,  that  we  did  not  test  with  the  20  pound  sledge  hammer 
at  either  site  #2  or  #3.  It  is  possible  that  in  our  desire  to  work  at  the 
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highest  possible  resolution,  we  prematurely  ignored  the  sledge 
hammer  at  sites  #2  and  #3. 


Since  the  data  are  of  relatively  high  quality  in  terms  of 
frequency  and  signal  to  noise  ratio,  but  reflections  were  not  readily 
visible,  it  is  our  conclusion  that  the  acoustic  impedance  contrast  at 
the  loess/conglomerate  interface  was  not  strong  enough  to  generate  a 
perceptible  reflection.  Three  different  locations  were  tested  with 
identical  results,  although  the  first  site  near  the  front  toe  of  the  dam 
had  data  contaminated  by  the  presence  of  the  slurry  wall.  The  lack 
of  acoustical  contrast  could  be  due  to  uneven  cementation  in  the 
conglomerate,  or  it  could  be  due  to  overall  bulk  velocity  being  near 
that  of  the  loess.  There  could  also  be  rough  surfaces  at  the  acoustical 
interface  that  inhibit  reflections,  although  no  diffraction  patterns 
were  observed  that  would  be  characteristic  of  some  varieties  of 
surface  roughness. 

Sufficient  information  is  available  to  calculate  the  approximate 
size  of  the  first  Fresnel  zone  which  gives  the  approximate  order  of 
direct  horizontal  resolution  that  would  be  possible  with  the  reflection 
method.  With  a  frequency  of  150  Hz,  a  velocity  of  3,000  feet/second, 
and  a  depth  of  100  feet,  the  diameter  of  the  Fresnel  zone  is  63  feet. 
For  the  same  velocity  and  frequency  at  a  depth  of  200  feet,  the 
Fresnel  zone  diameter  is  89  feet.  In  other  words,  direct  detection  of 
voids  less  than  perhaps  30  feet  across  at  this  site  at  depths  of  more 
than  100  feet  would  be  unlikely. 

The  figures  generated  for  this  report  are  in  a  walkaway  format. 
The  figure  descriptions  are  included  on  a  figure-by-figure  basis  and 
are  somewhat  redundant  when  compared  to  each  other.  The 
descriptions  will  however  allow  the  separation  of  individual  figures 
from  this  report  without  saciificing  the  basic  description  of  the 
acquisition  procedure  and  significant  information  about"  that  figure. 

No  digital  enhancement  techniques  were  performed  on  the 
seismic  data  other  than  variation  of  display  parameters. 
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SLEDGE  HAMMER  S1TE#1 
4  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  144  ft.  The  source  was  a  20 
lb  sledge  with  an  approximately  equivalent  weight  steel  plate.  The  4  hz 
lowcut  filters  were  the  lowest  low  cut  available  on  the  seismograph,  they 
therefore  represent  lowcuts  out.  The  geophone  spread  used  for  the  walkaway 
test  included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  In  order  to  obtain  source  receiver  offsets 
greater  than  96  ft  it  was  necessary  to  move  the  source  location.  Changing  the 
source  location  resulted  in  a  slight  trace  to  trace  mismatch  between  files 
recorded  with  different  source  locations.  This  can  be  seen  between  receivers 
with  offset  distances  of  46  ft  and  48  ft  on  this  72  trace  walkaway  plot.  The 
direct  and  refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  at  site#l, 
are  very  unusual  (in  amplitude,  trace  to  trace  phase,  and  relative  arrival 
pattern)  and  seem  to  indicated  significant  near-surface  structure.  The  direct 
wave  arrival  on  this  record  is  quite  substantial  with  respect  to  amplitude. 

The  wavelet  associated  with  the  direct  wave  appears  to  'ring'  throughout  the 
recording  period.  This  cyclic  nature  of  the  first  arrival  maybe  characteristic  of 
this  area.  The  direct  wave  velocity  at  this  site  appears  to  be  approximately 
1500  ft/sec.  The  apparent  first  refracting  velocity  is  approximately  4800  ft/sec. 
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SLEDGE  HAMMER  SITE*1 
4  HZ  ANALOG  LOWCUT  FILTER 
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SLEDGE  HAMMER  SITE#1 
128  HZ  ANALOG  LOWCUT  FILTER 

This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  144  ft.  The  source  was  a  20 
lb  sledge  with  an  approximately  equivalent  weight  steel  plate.  The  128  hz 
lowcut  filters  used  during  recording  of  this  walkaway  possess  a  18  dB/octave 
slope.  The  geophone  spread  used  for  this  walkaway  test  included  48  strings  of 
3  L-28E  40  Hz  geophones  Wired  in  series  and  planted  perpendicular  the  survey 
line.  In  order  to  obtain  source  receiver  offsets  greater  than  96  ft  it  was 
necessary  to  move  the  source  location.  Changing  the  source  location  resulted 
in  a  slight  trace  to  trace  mismatch  between  files  recorded  with  different  source 
locations.  This  can  be  seen  between  receivers  with  offset  distances  of  46  ft  and 
48  ft  on  this  72  trace  walkaway  plot.  The  direct  and  refracted  arrivals  on  this 
walkaway,  as  on  all  walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace 
to  trace  phase,  and  relative  arrival  pattern)  and  seem  to  indicated  significant 
near-surface  structure.  The  direct  wave  arrival  on  this  record  is  quite 
substantial  with  respect  to  amplitude.  The  wavelet  associated  with  the  direct 
wave  appears  to  'ring'  throughout  the  recording  period.  This  cyclic  nature  of 
the  first  arrival  maybe  characteristic  of  this  area.  The  direct  wave  velocity  at 
this  site  appears  to  be  approximately  1500  ft/sec.  The  apparent  first  refracting 
velocity  is  approximately  4800  ft/sec.  The  dominant  frequency  of  the  recorded 
signal  is  approximately  100  hz.  The  bandwidth  of  the  recorded  signal  is 
relatively  broad  with  frequencies  ranging  from  40  to  about  300  hz.  No  clearly 
identifiable  reflection  signal  is  present  on  this  walkaway. 
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SLEDGE  HAMMER  SITE^I 
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SLEDGE  HAMMER  SITE#1 
256  HZ  ANALOG  LOWCUT  FETER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  n:\aximum  source  to  receiver  offset  of  144  ft.  The  source  was  a  20 
lb  sledge  with  an  approximately  equivalent  weight  steel  plate.  The  256hz 
lowcut  filters  used  during  recording  of  this  walkaway  possess  a  18  dB/ octave 
slope.  The  geophone  spread  used  for  this  walkaway  test  included  48  strings  of 
3  L-28E  40  Hz  geophones  wired  in  series  and  planted  perpendicular  the  survey 
line.  In  order  to  obtain  source  receiver  offsets  greater  than  96  ft  it  was 
necessary  to  move  the  source  location.  Changing  the  source  location  resulted 
in  a  slight  trace  to  trace  mismatch  between  files  recorded  with  different  source 
locations.  This  can  be  seen  between  receivers  with  offset  distances  of  46  ft  and 
48  ft  on  this  72  trace  walkaway  plot.  The  direct  and  refracted  arrivals  on  this 
walkaway,  as  on  all  walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace 
to  trace  phase,  and  relative  arrival  pattern)  and  seem  to  indicated  significant 
near-surface  structure.  The  direct  wave  arrival  on  this  record  is  quite 
substantial  with  respect  to  amplitude.  The  wavelet  associated  with  the  direct 
wave  appears  to  'ring'  throughout  the  recording  period.  This  cyclic  nature  of 
the  first  arrival  maybe  characteristic  of  this  area.  Tne  direct  wave  velocity  at 
this  site  appears  to  be  approximately  1500  ft/sec.  The  apparent  first  refracting 
velocity  is  approximately  4800  ft/sec.  At  offsets  in  excess  of  100  ft  at 
approximately  100  msec  an  event  with  some  reflection  characteristics  can  be 
interpreted.  The  curvature  of  this  event  in  conjunction  with  its  arrival  time 
suggest  a  depth  of  approximately  170  ft.  The  air  wave  is  identifiable  on  this 
walkaway.  The  dominant  frequency  of  the  recorded  signal  is  approximately 
120  hz.  The  bandwidth  of  the  recorded  signal  is  relatively  broad  with 
frequencies  ranging  from  40  to  about  300  hz. 
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SLEDGE  HAMMER  S1TE#1 
384  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  144  ft.  The  source  was  a  20 
lb  sledge  with  an  approximately  equivalent  weight  steel  plate.  The  384  h2 
lowcut  filters  used  during  recording  of  this  walkaway  possess  a  18  dB/octave 
slope.  The  geophone  spread  used  for  this  walkaway  test  included  48  strings  of 
3  L-28E  40  Hz  geophones  wired  in  series  and  planted  perpendicular  the  survey 
line.  In  order  to  obtain  source  receiver  offsets  greater  than  96  ft  it  v;as 
necessary  to  move  the  source  location.  Changing  the  source  location  resulted 
in  a  slight  trace  to  trace  mismatch  between  files  recorded  with  different  source 
locations.  This  can  be  seen  between  receivers  with  offset  distances  of  46  ft  and 
48  ft  on  this  72  trace  walkaway  plot.  The  direct  and  refracted  arrivals  on  this 
walkaway,  as  on  all  walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace 
to  trace  phase,  and  relative  arrival  pattern)  and  seem  to  indicated  significant 
near-surface  structure.  The  direct  wave  arrival  on  this  record  is  quite 
substantial  with  respect  to  amplitude.  The  wavelet  associated  with  the  direct 
wave  appears  to  'ring*  throughout  the  recording  period.  This  cyclic  nature  of 
the  first  arrival  maybe  characteristic  of  this  area.  The  direct  wave  velocity  at 
this  site  appears  to  be  approximately  1500  ft/sec.  The  apparent  first  refracting 
velocity  is  approximately  4800  ft/sec.  At  offsets  in  excess  of  100  ft  at 
approximately  100  msec  an  event  with  some  reflection  characteristics  can  be 
interpreted.  The  curvature  of  this  event  in  conjunction  with  its  arrival  time 
suggest  a  depth  of  approximately  170  ft.  The  air  wave  is  identifiable  on  this 
walkaway.  The  dominant  frequency  of  the  recorded  signal  is  approximately 
150  hz.  The  bandwidth  of  the  recorded  signal  is  relatively  broad  with 
frequencies  ranging  from  40  to  about  300  hz.  The  high  frequency  (600  hz) 
horizontally  coherent  spikes  obvious  above  90  msec  at  offsets  greater  than  90 
ft  area  equipment  related  and  do  not  represent  recorded  acoustical  energy. 
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SLEDGE  HAMMER  SITE#1 
512  HZ  ANAl-OG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  with  minimum  source  to  receiver  offset  of  48  ft  and  maximum  source 
to  receiver  offset  of  144  ft.  The  source  was  a  20  lb  sledge  with  an 
approximately  equivalent  weight  steel  plate.  The  512  hz  lowcut  filters  used 
during  recording  of  this  walkaway  possess  a  18  dB/octave  slope.  The 
geophone  spread  used  for  this  walkaway  test  included  48  strings  of  3  L-28E  40 
Hz  geophones  wired  in  series  and  planted  perpendicular  the  survey  line. 

The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  at 
site#l,  are  very  unusual  (in  amplitude,  trace  to  trace  phase,  and  relative 
arrival  pattern)  and  seem  to  indicated  significant  near-surface  structure.  The 
direct  wave  arrival  on  this  record  is  quite  substantial  w'ith  respect  to 
amplitude.  The  wavelet  associated  with  the  direct  wave  appears  to  'ring' 
throughout  the  recording  period.  This  cyclic  nature  of  the  first  arrival  maybe 
characteristic  of  this  area.  The  direct  wave  velocity  at  this  site  appears  to  be 
approximately  1500  ft /sec.  The  apparent  first  refracting  velocity  is 
approximately  4800  ft/sec.  At  offsets  in  excess  of  100  ft  at  approximately  100 
msec  an  event  with  some  reflection  characteristics  can  be  interpreted.  The 
curvature  of  this  event  in  conjunction  with  its  arrival  time  suggest  a  depth  of 
approximately  170  ft.  The  air  wave  is  identifiable  on  this  walkaway.  The 
dominant  frequency  of  the  recorded  signal  is  approximately  150  hz.  The 
bandwidth  of  the  recorded  signal  is  relatively  broad  with  frequencies  ranging 
from  40  to  about  300  hz.  The  increased  amounts  of  high  frequency  noise 
present  on  the  early  portions  of  the  record,  results  from  the  removal  of  a 
significant  amount  of  lower  frequency  signal  by  the  severe  lowcut  filter.  This 
severe  filtering  allows  low  amplitude  noise  to  become  a  more  significant 
portion  of  the  data  set. 
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DOWNHOLE  30.06  SHE#! 

4  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  144  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  4  hz  lowcut  filter.s  were 
the  lowest  low  cut  available  on  the  seismograph,  they  therefore  represent 
lowcuts  out.  The  geophone  spread  used  for  the  walkaway  test  included  48 
strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted  perpendicular 
the  survey  line.  In  order  to  obtain  source  receiver  offsets  greater  than  96  ft  it 
was  necessary  to  move  the  source  location.  Changing  the  source  location 
resulted  in  a  slight  trace  to  trace  mismatch  between  files  recorded  with 
different  source  locations.  This  can  be  seen  between  receivers  with  offset 
distances  of  94  ft  and  96  ft  on  this  72  trace  walkaway  plot.  The  direct  and 
refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  at  site#l,  are  very 
unusual  (in  amplitude,  trace  to  trace  phase,  and  relative  arrival  pattern)  and 
seem  to  indicated  significant  near-surface  structure.  The  direct  wave  arrival 
on  this  record  is  quite  substantial  with  respect  to  amplitude.  The  wavelet 
associated  with  the  direct  wave  appears  to  'ring'  throughout  the  recording 
period.  This  cyclic  nature  of  the  first  arrival  maybe  characteristic  of  this  area. 
The  direct  wave  velocity  at  this  site  appears  to  be  approximately  1500  ft/sec. 
The  apparent  first  refracting  velocity  is  approximately  4800  ft/sec.  The  4  hz 
walkaway  with  the  sledge  hammer  and  steel  plate  is  quite  consistent  with  this 
4  hz  walkaway.  The  main  distinction  is  the  frequency  of  the  identifiable  body 
wave  energy.  The  30.06  appears  to  have  a  dominant  frequency  at  least  30  to  40 
hz  higher. 


A17 


TIME  (MSEC) 


DOWNHOLE  30.06  SITE^l 
4  HZ  ANALOG  LOWCUT  FILTER 


SOURCE  TO  RECEIVER  OFFSET 

(FEET) 

0  144 


0  48 


I _ I 

DISTANCE 

(FEET) 


A18 


DOWNHOLE  30.06  SITE#! 

128  HZ  ANALOG  LOWCUT  HLTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  128  hz  lowait  filters 
provide  a  18  dB/ octave  rolloff.  TTie  geophone  spread  used  for  the  walkaway 
test  included  48  strings  of  3  L‘‘28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  The  apparent  phase  anomaly  between 
receivers  at  46  ft  and  48  ft  are  the  result  of  source  static.  The  30.06  rifle  xv  as 
fired  into  the  same  hole  for  several  shots.  This  procedure  in  some 
environments  can  result  in  a  slight  time  discrepancy  on  multi'file 
walkaways.  The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all 
walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace  to  trace  phase,  and 
relative  arrival  pattern)  and  seem  to  indicated  significant  near-surface 
structure.  The  direct  wave  arrival  on  this  record  is  quite  substantial  with 
respect  to  amplitude.  The  wavelet  associated  with  the  direct  wave  appears  to 
'ring'  throughout  the  recording  period.  This  cyclic  nature  of  the  first  arrival 
maybe  characteristic  of  this  area.  The  direct  wave  velocity  at  this  site  appears 
to  ^  approximately  1500  ft/sec.  The  apparent  first  refracting  velocity  is 
approximately  4800  ft/sec.  No  obvious  primary  reflection  energy  can  be 
interpreted  on  this  walkaway. 
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DOWNHOLE  30.06  SITE#1 
256  HZ  ANALOG  LOWCUT  FILTER 


This  filie  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  256  hz  lowcut  filters 
provide  a  18  dB/octave  rolloff.  The  geophone  spread  used  for  the  walkaway 
test  included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  The  apparent  phase  anomaly  between 
receivers  at  46  ft  and  48  ft  are  the  result  of  source  static.  The  30.06  rifle  was 
fired  into  the  same  hole  for  several  shots.  This  procedure  in  some 
environments  can  result  in  a  slight  time  discrepancy  on  multi-file 
walkaways.  The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all 
walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace  to  trace  phase,  and 
relative  arrival  pattern)  and  seem  to  indicated  significant  near-surface 
structure.  The  direct  wave  arrival  on  this  record  is  quite  substantial  with 
respect  to  amplitude.  The  wavelet  associated  with  the  direct  wave  appears  to 
'ring'  throughout  the  recording  period.  This  cyclic  nature  of  the  first  arrival 
maybe  characteristic  of  this  area.  The  direct  wave  velocity  at  this  site  appears 
to  be  approximately  1500  ft/sec.  The  apparent  first  refracting  velocity  is 
approximately  4800  ft/sec.  There  appears  to  be  more  coherent  energy  above 
the  airwave  and  below  the  refracted  arrival  on  the  30.06  data  in  comparison  to 
equivalent  sledge  hammer  data.  The  frequency  of  the  30.06  body  wave 
information  possess  a  slightly  higher  frequency.  No  obvious  primary 
reflection  energy  can  be  interpreted  on  this  walkaway. 
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DOWNHOLE  30.06  SITE#! 

512  HZ  ANALOG  LOWCLTT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  512  hz  lowcut  filters 
provide  a  18  dB/octave  rolloff.  TTje  geophone  spread  used  for  the  walkaway 
test  included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  The  apparent  phase  anomaly  between 
receivers  at  46  ft  and  48  ft  are  the  result  of  source  static.  The  30.06  rifle  was 
fired  into  the  same  hole  for  several  shots.  This  procedure  in  some 
environments  can  result  in  a  slight  time  discrepancy  on  multi-file 
walkaways.  The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all 
walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace  to  trace  phase,  and 
relative  arrival  pattern)  and  seem  to  indicated  significant  near-surface 
structure.  The  direct  wa^'e  arrival  on  this  record  is  quite  substantial  with 
respect  to  amplitude.  The  wavelet  associated  with  the  direct  wave  appears  to 
'ring*  throughout  the  recording  period.  This  cyclic  nature  of  the  first  arrival 
maybe  characteristic  of  this  area.  The  direct  wave  velocity  at  this  site  appears 
to  be  approximately  1500  ft/sec.  The  apparent  first  refracting  velocity  is 
approximately  4800  ft/sec.  There  appears  to  be  more  coherent  energy  above 
the  airwave  and  below  the  refracted  arrival  on  the  30.06  data  in  comparison  to 
equivalent  sledge  hammer  data.  The  frequency  of  the  30.06  body  wave 
information  possess  a  slightly  higher  frequency.  There  appears  to  be  little 
difference  between  the  256  lowcut  and  the  512  lowcut  except  for  a  slightly 
higher  amount  of  electronic  noise  evident  pre-first  arrival  on  the  more 
distant  set  of  24  channels  (48-96  ft).  No  obvious  primary  reflection  energy  can 
be  interpreted  on  this  walkaway. 
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DOWNHOLE  so  CAL  SITE#1 
384  HZ  ANALOG  LOWCUT  FILTER 


Thus  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  50  cal  rifle  firing  a  750  grain  projectile.  The  384  hz  lowcut  filters 
provide  a  18  dB/ octave  rolloff.  The  geophone  spread  used  for  the  walkaway 
test  included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  The  apparent  phase  anomaly  between 
receivers  at  46  ft  and  48  ft  are  the  result  of  source  static.  The  50  cal  rifle  was 
fired  into  the  same  hole  for  several  shots.  This  procedure  in  some 
environments  can  result  in  a  slight  time  discrepancy  on  multi-file 
walkaways.  The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all 
walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace  to  trace  phase,  and 
relative  arrival  pattern)  and  seem  to  indicated  significant  near-surface 
structure.  The  direct  wave  arrival  on  this  record  is  quite  substantial  with 
respect  to  amplitude.  The  wavelet  associated  with  the  direct  wave  appears  to 
'ring'  throughout  the  recording  period.  This  cyclic  nature  of  the  first  arrival 
maybe  characteristic  of  this  area.  The  direct  wave  velocity  at  this  site  appears 
to  be  approximately  1500  ft/sec.  The  apparent  first  refracting  velocity  is 
approximately  4800  ft/sec.  There  appears  to  be  more  coherent  energy  above 
the  airwave  and  below  the  refracted  arrival  on  the  50  cal  data  in  comparison 
to  either  equivalent  sledge  hammer  or  30.06  data.  A  significantly  larger 
amount  of  total  energy  is  generated  by  the  SO  cal  than  either  of  the  other  2 
source  tried  at  this  location.  No  obvious  primary  reflection  energy  can  be 
interpreted  on  this  walkaway. 
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DOWNHOLE  so  CAL  SITE#! 

512  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  50  cal  rifle  firing  a  750  grain  projectile.  The  512  hz  lowcut  filters 
provide  a  18  dB/octave  rolloff.  The  geophone  spread  used  for  the  walkaway 
test  included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  The  apparent  phase  anomaly  betv\>^een 
receivers  at  46  ft  and  48  ft  are  the  result  of  source  static.  The  50  cal  rifle  was 
fired  into  the  same  hole  for  several  shots.  This  procedure  in  some 
environments  can  result  in  a  slight  time  discrepancy  on  multi-file 
walkaways.  The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all 
walkaways  at  site#l,  are  very  unusual  (in  amplitude,  trace  to  trace  phase,  and 
relative  arrival  pattern)  and  seem  to  indicated  significant  near-surface 
structure.  The  direct  wave  arrival  on  this  record  is  quite  substantial  with 
respect  to  amplitude.  The  wavelet  associated  with  the  direct  wave  appears  to 
'ring*  throughout  the  recording  period.  This  cyclic  nature  of  the  first  arrival 
maybe  characteristic  of  this  area.  The  direct  wave  velocity  at  tliis  site  appears 
to  be  approximately  1500  ft/sec.  The  apparent  first  refracting  velocity  is 
approximately  4800  ft/sec.  There  appears  to  be  more  coherent  energy  above 
the  airwave  and  below  the  refracted  arrival  on  the  50  cal  data  in  comparison 
to  either  equivalent  sledge  hammer  or  30.06  data.  A  significantly  larger 
amount  of  total  energy  is  generated  by  the  50  cal  than  either  of  the  other  2 
source  tried  at  this  location.  Most  of  the  arrivals  within  the  time  window 
between  the  refractions  and  the  ground  roll  have  a  slope  very  similar  to 
either  the  refractions  or  the  direct  wave.  Making  the  likelihood  of  any  of 
those  events  being  reflections  relatively  poor.  No  obvious  primary  reflection 
energy  can  be  interpreted  on  this  walkaway. 


A27 


DOWNHOLE  30.06  SITE#! 

4  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  144  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  4  hz  lowcut  filters  were 
the  lowest  low  cut  available  on  the  seismograph,  they  therefore  represent 
lowcuts  out.  The  geophone  spread  used  for  the  walkaway  test  included  48 
strings  of  3  L--28E  40  i\z  geophones  wired  in  series  and  planted  perpendicular 
the  survey  line.  In  order  to  obtain  source  receiver  offsets  greater  than  96  ft  it 
was  necessary  to  move  the  source  location.  Changing  the  source  location 
resulted  in  a  slight  trace  to  trace  mismatch  between  files  recorded  with 
different  source  locations.  This  can  be  seen  between  receivers  with  offset 
distances  of  94  ft  and  96  ft  on  this  72  trace  walkaway  plot.  The  direct  and 
refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  from  this  site  and  site 
#3  are  much  more  characteristic  of  “normal''  first  arrivals.  The  direct  wave 
velocity  at  this  site  appears  to  be  approximately  750  ft/sec.  The  apparent  first 
refracting  velocity  is  approximately  1400  ft/sec.  As  at  site#l,  the  majority  of 
the  energy  arriving  within  the  classically  quite  window  between  the 
refraction  and  ground  roll  arrivals  have  a  slope  very  similar  to  that  of  first 
arrival  refraction.  As  well,  the  NMO  velocity  of  these  arrivals,  clearly 
disqualify  them  from  consideration  as  reflection  energy.  No  reflection 
events  can  be  confidently  identified  on  this  seismogram. 
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DOWNHOLE  30.06  SITE#2 
128  HZ  ANALOG  LOWCUT  FILTER 

This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  128  hz  lowcut  filters 
have  an  18  dB /octave  slope.  The  geophone  spread  used  for  the  walkaway  test 
included  48  strings  of  3  L~28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  A  slight  mismatdt  can  be  observed  between 
offsets  46  and  48  as  a  result  of  different  recorded  shots.  The  same  hole  was 
occupied  for  both  shot  but  the  second  shot  in  the  hole  will  have  a  slightly 
different  travel  path  than  the  first.  This  different  travel  path  is  mainly 
represented  as  deeper  penetration.  The  direct  and  refracted  arrivals  on  this 
walkaway,  as  on  all  walkaways  from  this  site  and  site  #3  are  much  more 
characteristic  of  "normal"  first  arrivals.  The  direct  wave  velocity  at  this  site 
appears  to  be  approximately  750  ft/ sec.  The  apparent  first  refracting  velocity  is 
approximately  1400  ft/sec.  As  at  site#l,  the  majority  of  the  energy  arriving 
within  the  classically  quite  window  between  the  refraction  and  ground  roll 
arrivals  have  a  slope  very  similar  to  that  of  first  arrival  refraction.  As  well, 
the  NMO  velocity  of  these  arrivals,  clearly  disqualify  them  from 
consideration  as  reflection  energy.  No  reflection  events  can  be  confidently 
identified  on  this  seismogram.  The  ground  roll  arrivals  at  offsets  less  than  15 
ft  lack  coherency.  This  is  probably  related  to  the  low  cut  filtering  of  the 
systems  filters  and  the  40  hz  geophones. 
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DOWNHOLE  30.06  SITE#2 
256  HZ  ANALOG  LOWCUT  HLTER 

This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  256  hx  lowcut  filters 
have  an  18  dB /octave  slope.  The  geophone  spread  used  for  the  walkaway  test 
included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  A  slight  mismatch  can  be  observed  between 
offsets  46  and  48  as  a  result  of  different  recorded  shots.  The  same  hole  was 
occupied  for  both  shot  but  the  second  shot  in  the  hole  will  have  a  slightly 
different  travel  path  than  the  first.  This  different  travel  path  is  mainly  the 
result  of  deeper  penetration  of  the  second  shot  in  comparison  to  the  first. 

The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  from 
this  site  and  site  #3  are  much  more  characteristic  of  "normal"  first  arrivals. 
The  direct  wave  velocity  at  this  site  appears  to  be  approximately  750  ft/sec. 

The  apparent  first  refracting  velocity  is  approximately  1400  ft/sec.  As  at 
site#l,  the  majority  of  the  energy  arriving  within  the  classically  quite  window 
between  the  refraction  and  ground  roll  arrivals  have  a  slope  very  similar  to 
that  of  first  arrival  refraction.  As  well,  the  NMO  velocity  of  these  arrivals, 
clearly  disqualify  them  from  consideration  as  reHection  energy.  No  reflection 
events  can  be  confidently  identified  on  this  seismograni.  The  ground  roll 
arrivals  at  offsets  less  than  15  ft  lack  coherency.  This  is  probably  related  to  the 
low  cut  filtering  of  the  systems  filters  and  the  40  hz  geophones. 
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DOWNHOLE  30.06  SrrE#2 
384  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  384  hz  lowcut  filters 
have  an  18  dB/octave  slope.  The  geophone  spread  used  for  the  walkaway  test 
included  48  strings  of  3  L-2dE  40  Hz  geophones  wired  in  series  and  plant^ 
perpendicular  the  survey  line.  A  slight  mismatch  can  be  observed  between 
offsets  46  and  48  as  a  result  of  different  recorded  shots.  The  same  hole  was 
occupied  for  both  shot  but  the  second  shot  in  the  hole  will  have  a  slightly 
different  travel  path  than  the  first.  This  different  travel  path  is  mainly  the 
result  of  deeper  penetration  of  the  second  shot  in  comparison  to  the  first. 

The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  from 
this  site  and  site  #3  are  much  more  characteristic  of  "normal"  first  arrivals. 
The  direct  wave  velocity  at  this  site  appears  to  be  approximately  750  ft/sec. 

The  apparent  first  refracting  velocity  is  approximately  1400  ft/sec.  As  at 
site#l,  the  majority  of  the  energy  arriving  within  the  classically  quite  window 
between  the  refraction  and  ground  roll  arrivals  have  a  slope  very  similar  to 
that  of  first  arrival  refraction.  As  well,  the  NMO  velocity  of  these  arrivals, 
clearly  disqualify  them  from  consideration  as  reflection  energy.  The 
dominant  frequency  of  the  body  wave  energy  is  noticeably  higher  on  this 
record  than  on  the  previous  records  with  lower  lowcut  filters  from  this 
location.  No  reflection  events  can  be  confidently  identified  on  this 
seismogram.  The  ground  roll  arrivals  at  offsets  less  than  15  ft  lack  coherency. 
This  is  probably  related  to  the  low  cut  filtering  of  the  systems  filters  and  the  40 
hz  geophones. 
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DOWNHOLE  30.06  SrrE#2 
512  HZ  ANALOG  LOWCUT  FILTER 

This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft,  The  source  was  a 
downhole  50  cal  rifle  firing  a  750  grain  projectile.  The  512  hz  lowcut  filters 
have  an  18  dB/octave  slope.  The  geophone  spread  used  for  the  walkaway  test 
included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  A  slight  mismatch  can  be  observed  between 
offsets  46  and  48  as  a  result  of  different  recorded  shots.  The  same  hole  was 
occupied  for  both  shots  that  make  up  this  walkaway.  The  second  shot 
(recording  the  longer  offsets)  in  the  hole  had  slightly  deeper  penetration  in 
comparison  to  the  first.  The  direct  and  refracted  arrivals  on  this  walkaway,  as 
on  all  walkaways  from  this  site  and  site  #3  are  much  more  characteristic  of 
"normal"  first  arrivals  patterns.  The  direct  wave  velocity  at  this  site  appears 
to  be  approximately  750  ft/sec.  The  apparent  first  refracting  velocity  is 
approximately  1400  ft/sec.  As  at  site#l,  the  majority  of  the  energy  arriving 
within  the  classically  quite  window  between  the  refraction  and  ground  roll 
arrivals  have  a  slope  very  similar  to  that  of  first  arrival  refraction.  As  well, 
the  NMO  velocity  of  these  arrivals,  clearly  disqualify  them  from 
consideration  as  reflection  energy.  No  reflection  events  can  be  confidently 
identified  on  this  seismogram.  The  50  cal  clearly  generates  more  recordable 
seismic  energy  than  either  the  30.06  or  the  sledge  hammer  and  plate.  As  with 
some  of  the  earlier  walkaways,  pre-first  arrival  noise  can  be  interpreted.  The 
frequency  and  coherency  of  the  energy  is  indicative  of  instrument  noise.  The 
ground  roll  arrivals  at  offsets  less  than  15  ft  lack  trace  to  trace  coherency.  This 
is  probably  related  to  the  low  cut  filtering  of  the  systems  filters  and  the  40 
geophones.  Coherent  energy  arrivals  beyond  48  ft  of  offset  and  deeper  in  time 
than  approximately  100  msec  on  the  seismogram  possess  a  variety  of  phase 
velocities. 
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DOWNHOLE  30.06  SrrE#3 
384HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  144  ft.  The  source  was  a 
downhole  30,06  rifle  firing  a  180  grain  projectile.  The  384  hz  lowcut  filters 
have  an  18  dB/octave  slope.  The  geophone  spread  used  for  the  walkaway  test 
included  48  strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  A  slight  mismatch  can  be  observed  between 
offsets  46  and  48  as  a  result  of  different  recorded  shots.  The  same  hole  was 
occupied  for  the  first  and  second  shots.  The  shot  recorded  at  the  longer  offset 
had  a  slightly  different  depth  of  penetration  than  the  first  shot  in  the  hole. 

The  slight  phase  mismatch  between  the  traces  at  96  and  98  ft  of  offset  is  a 
result  of  the  different  source  location  and  therefore  a  different  source  static. 

The  direct  and  refracted  arrivals  on  this  walkaway/  as  on  all  walkaways  from 
this  site  are  much  more  characteristic  of  "normal"  first  arrivals.  The  direct 
wave  velocity  at  this  site  appears  to  be  approximately  625  ft/sec.  The  apparent 
first  refracting  velocity  is  approximately  1900  ft/sec.  As  at  site#l,  the  majority 
of  the  energy  arriving  within  the  classically  quite  window  between  the 
refraction  and  ground  roll  arrivals  have  a  slope  very  similar  to  that  of  first 
arrival  refraction.  As  well/  the  NMO  velocity  of  these  arrivals/  dearly 
disqualify  them  from  consideration  as  reflection  energy.  No  reflection 
events  can  be  confidently  identified  on  this  seismogram.  The  high  frequency 
spikes  on  the  inside  24  traces  within  the  first  100  msec  is  the  result  of 
overdriving  the  analog  filters.  The  operator  of  the  seismograph  has  no  way 
of  knowing  if  the  filters  are  overdriven  except  through  visual  inspection  of 
field  plots.  The  event  previously  identified  at  190  msec  as  a  potential  reflected 
event  is  present  but  much  more  obscured  on  this  walkaway  than  of  the  256  hz 
lowcut  display.  The  probably  of  this  event  being  a  reflection  is  neither 
enhanced  or  reduced  as  a  result  of  analysis  of  this  walkaway.  The  ground  ro 
arrivals  at  offsets  less  than  15  ft  lack  coherency.  This  is  probably  related  to  the 
low  cut  filtering  of  the  systems  filters  and  the  40  hz  geophones.  Three 
distinctly  different  phase  velodties  can  be  identified. 
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DOWNHOLE  30.06  SrrE#3 
256  HZ  ANALOG  LOWCUT  HLTER 

This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  96  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  256  hz  lowcut  filters 
have  an  18  dB/octave  slope.  T1\e  geophone  spread  used  for  the  walkaway  test 
included  48  strings  of  3  L‘28E  40  Hz  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  A  slight  mismatch  can  be  observed  between 
offsets  46  and  48  as  a  result  of  different  recorded  shots.  The  same  hole  was 
occupied  for  both  shots.  The  second  shot  in  the  hole  penetrated  deeper,  as  a 
direct  result  of  the  first  shots,  intuitively  obvious  deepening  of  the  hole.  The 
deeper  hole  allowed  the  second  projectile  to  dissipate  it's  energy  into  material 
possessing  a  slightly  different  velocity  as  well  as  travel  longer  in  the  air  of  the 
borehole.  The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all 
walkaways  from  this  site  are  much  more  characteristic  of  "normal"  first 
arrivals.  The  direct  wave  velocity  at  this  site  appears  to  be  approximately  600 
ft/sec.  The  apparent  first  refracting  velocity  is  approximately  1800  ft/sec.  As 
at  site#l,  the  majority  of  the  energy  arriving  within  the  classically  quite 
window  between  the  refraction  and  ground  roll  arrivals  have  a  slope  very 
similar  to  that  of  first  arrival  refraction.  As  well,  the  NMO  velocity  of  these 
arrivals,  clearly  disqualify  them  from  consideration  as  reflection  energy.  A 
subtle  hint  of  coherency  can  be  identified  at  approximately  190  msec  on  th 
inside  35  traces.  This  event  possess  to  high  a  velocity  at  these  offsets  and  trace 
spacings  to  confidently  calculate  the  NMO  velocity.  With  some  spectral 
analysis  and  digital  filtering  this  event  could  possibly  be  enhanced.  Some 
possibility  of  spatial  aliasing  exists  with  this  event.  *1116  only  way  to 
confidently  eliminate  that  possibility  is  by  varying  the  source  offset  slightly. 
This  of  course  is  not  an  option  at  this  time.  An  F-k  filter  probably  would 
enhance  the  coherency  of  this  event.  The  frequency  and  lack  of  good 
coherency  suggests  caution  be  taken  in  confidently  identifying  this  event  as  a 
reflection. 
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DOWNHOLE  30.06  SITE#3 
128  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maxinaum  source  to  receiver  offset  of  144  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  128  hz  lowcut  filters 
have  an  18  dB/octave  slope.  The  geophone  spread  used  for  the  walkaway  test 
included  48  strings  of  3  L-28E  40  geophones  wired  in  series  and  planted 
perpendicular  the  survey  line.  A  slight  mismatch  can  be  observed  between 
offsets  46  and  48  as  a  result  of  different  recorded  shots.  The  same  hole  was 
occupied  for  the  first  and  second  shots.  The  shot  recorded  at  the  longer  offset 
had  a  slightly  different  depth  of  penetration  than  the  first  shot  in  the  hole. 
The  slight  phase  mismatch  between  the  traces  at  96  and  98  ft  of  offset  is  a 
result  of  the  different  source  location  and  therefore  a  different  source  static. 
The  direct  and  refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  from 
this  site  are  much  more  characteristic  of  "normal"  first  arrivals.  The  direct 
wave  velocity  at  this  site  appears  to  be  approximately  625  ft/sec.  The  apparent 
first  refracting  velocity  is  approximately  1900  ft/sec.  As  at  site#l,  the  majority 
of  the  energy  arriving  within  the  classically  quite  window  between  the 
refraction  and  ground  roll  arrivals  have  a  slope  very  similar  to  that  of  first 
arrival  refraction.  As  well,  the  NMO  velocity  of  these  arrivals,  clearly 
disqualify  them  from  consideration  as  reflection  energy.  No  reflection 
events  can  be  confidently  identified  on  this  seismogram.  The  ground  roll 
arrivals  at  offsets  less  than  15  ft  lack  coherency.  This  is  probably  related  to  the 
low  cut  filtering  of  the  systems  filters  and  the  40  hz  geophones.  Three 
distinctly  different  phase  velocities  can  be  identified. 
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DOWNHOLE  30.06  SrrE#3 
4  HZ  ANALOG  LOWCUT  FILTER 


This  file  represents  a  walkaway  noise  test  with  receivers  on  2  foot 
centers  and  maximum  source  to  receiver  offset  of  144  ft.  The  source  was  a 
downhole  30.06  rifle  firing  a  180  grain  projectile.  The  4  hz  lowcut  filters  were 
the  lowest  low  cut  available  on  the  seismograph;  they  therefore  represent 
lowcuts  out.  The  geophone  spread  used  for  the  walkaway  test  included  48 
strings  of  3  L-28E  40  Hz  geophones  wired  in  series  and  planted  perpendicular 
the  survey  line.  In  order  to  obtain  source  receiver  offsets  greater  than  96  ft  it 
was  necessary  to  move  the  source  location.  Changing  the  source  location 
resulted  in  a  slight  trace  to  trace  mismatch  between  files  recorded  with 
different  source  locations.  This  can  he  seen  between  receivers  with  offset 
distances  of  94  ft  and  96  ft  on  this  72  trace  walkaway  plot.  The  direct  and 
refracted  arrivals  on  this  walkaway,  as  on  all  walkaways  from  this  site  are 
much  more  characteristic  of  "normal"  first  arrivals.  The  direct  wave  velocity 
at  this  site  appears  to  be  approximately  650  ft/sec.  The  apparent  first  refracting 
velocity  is  approximat.:.iy  1800  ft/sec.  As  at  site#l,  the  majority  of  the  energy 
arriving  within  the  classically  quite  window  betw'een  the  refraction  and 
ground  roll  arrivals  have  a  slope  very  similar  to  that  of  first  arrival  refraction. 
As  well,  the  NMO  velocity  of  these  arrivals,  dearly  disqualify  them  from 
consideration  as  reflection  energy.  Tlie  energy  arriving  in  the  window  most 
likely  to  possess  Hrst  order  reflection  ^formation  is  relatively  broad  band. 

The  frequency  band  extends  from  about  40  hz  through  350  hz.  No  reflection 
events  can  be  confidently  identified  on  this  seismogram. 
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APPENDIX  B 

GROUND  PENETRATING  RADAR  SURVEY,  MILL  CREEK  DAM 

Wllllanson  and  Associates 
Seattle,  Washington 


RESULTS  OP  A  GPR  SURVEY 
AT  TBE 

If  ILL  CREEK  DAM.  WALLA  WALLA.  WASHINGTON 


INTRODUCTION 

This  report  presents  the  results  of  a  geophysical  survey,  conducted  for  the 
USACOB  at  Mill  Creek  Dam,  Walla  Walla,  Washington.  The  field  work,  using 
ground  penetrating  radar  (GPR),  was  performed  on  May  2nd  and  3rd,  1989. 

The  objective  of  the  survey  was  to  detect  cavities  in  a  sUt  deposit  which  is 
the  foundation  of  the  dam  as  well  as  the  right  abutment.  The  silt  is  reported 
to  be  from  40  to  100  feet  thick.  Beneath  the  silt  is  a  conglomerate  that  is 
approximately  100  feet  thick  which  rests  on  basalt.  Sinkholes  have  devel-* 
oped  in  the  silts  on  the  floor  of  the  reservoir  and  there  is  concern  that  cav¬ 
ities  in  the  silt  might  threaten  the  integrity  of  the  dam.  A  suggested  mech¬ 
anism  for  cavity  formation  is  by  piping  d*  sediments  down  through  the 
rather  pervious  zones  within  the  conglomerate. 

GPR  TECHNIQUE 

The  GPR  produces  a  continuous  prdile  of  the  subsurface  soils  and  bedrock 
features  along  the  survey  traverse.  This  is  accomplished  by  towing  a  radar 
antenna,  at  a  walking  pace,  along  the  desired  trackline.  The  antenna  is  con¬ 
nected  by  a  40  meter  cable  to  a  processor  and  display  recorder  mounted  in 
the  field  vehicle  (Figure  1 ).  During  the  traverse,  the  antenna  alternately 
transmits  and  receives  high  frequency  (120  MHz)  radar  pulses.  The  received 
pulses,  which  represent  reflections  from  subsurface  features  are  immediate¬ 
ly  processed  and  displayed  on  the  graphic  recorder.  Since  this  transmit/rec¬ 
eive  process  occurs  at  such  a  high  rate,  a  vertical  profile  or  "sounding "  is 
obtained  for  each  4  inches  of  horizontal  travel  along  the  survey  transect. 

The  resulting  display  is  a  continuous  subsurface  profile  providing  detailed 
information  on  subsurface  features.  Those  features  on  the  radargram  that 
represent  soil  horizons  or  stratigraphic  layers  are  produced  by  changes  in 
the  conductivity  between  adjacent  layers  of  soil.  Typically  these  variations 
can  be  associated  with  changes  in  lithology,  an  increase  in  moisture  content, 
or  contaminants  having  electrical  conductivity  properties  different  than 
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Figure  i  Illustration  of  GPR  technique  and  radargram. 
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those  of  water.  lo  addition,  reflections  can  be  produced  by:  subtle  variations 
in  the  geologic  structure  and  fabric  such  as  fractures,  faulting  and  subsurface 
voids. 

FIELD  PROCEDURES 

The  survey  was  performed  with  a  GSSI  System  3  ground  penetrating  radar 
(GPR)  system  using  a  120  MH2  antenna.  Data  was  collected  along  four  pri¬ 
mary  survey  lines,  selected  by  the  ACOE,  and  around  an  eiisting  sinkhole  in 
the  bottom  of  the  reservoir,  ^veral  runs  were  made  along  each  line  with 
the  display  recorder  set  for  different  scales  on  each  run  in  an  attempt  to 
maximize  the  depth  of  penetration  or  the  near  surface  resolution.  A  1 50 
nanosecond  scale  was  used  to  provide  detailed  information  on  the  near 
surface  stratigraphy  in  the  area  of  the  sinkhole  while  a  300  nanosecond  scale 
was  used  to  map  deeper  subsurface  features  on  lines  C.  D.  E  &  F.  These  time 
scales  represent  a  depth  scale  of  approximately  30  feet  and  60  feet. 

The  antenna  was  pulled  by  hand  along  each  of  the  survey  lines  that  had 
been  previously  laid  out  and  staked  by  the  ACOE.  As  the  antenna  was  pul¬ 
led  past  each  stake  (20  foot  intervals)  a  control  or  position  mark  was  placed 
on  the  radargram.  Pinal  position  of  the  lines  and  horizontal  control  was  pro¬ 
vided  by  the  ACXIE  for  the  final  report. 

Field  notes  describing  the  general  condition  of  the  ground  surface,  such  as 
major  topographic  changes,  the  presence  of  cultural  features  or  waste  mat¬ 
erials,  and  the  type  of  surficial  soils  were  kept  during  the  survey.  This  in¬ 
formation  used  was  to  assist  in  analyzing  the  GPR  data. 

DATA  ANATSIS  AND  INTERPRETATION; 

The  areas  surveyed  with  the  GPR  are  plotted  on  Plate  1  which  shows  both 
the  tracklines  and  the  location  of  subsurface  anomalies.  The  anomaly  num¬ 
bers  on  this  plate  correspond  with  the  numbers  shown  on  the  radargrams 
located  in  the  appendix.  A  brief  discussion  of  the  subsurface  geology  and 
significant  features  is  presented  in  the  following  table. 
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LOfflaiy^aber.  Description 

General  Thin  surface  layer  ( 1  to  3  feet  thick)  along  entire 

survey  line.  Below  this  reflector  materials  are 
homogeneous  with  only  an  occassional  subsurface 
reflector. 


C-42  Subsurface  ringing  from  -2  feet  to  -1 0  feel 

C-^9  to  C-63  Shallow  subsurface  depression 

C-60  Smalt  diffractions  at  - 17  feet,  some  shallower  re¬ 

flectors. 


LINED 


tomaly  Number 


DgRCtiptioa 


General 
D-51  to  D-53 


Thin  surface  layer  (1  to  3  feet  thick)  and  very 
homogenous  and  uniform  subsurface  materials. 

Small  subsurface  depression 


LINEB 

LamaiyUumbef,  Description 

General  Thin  surface  layer  ( 1  to  3  feet  thick)  and  homo¬ 

genous  subsurface  sediments. 


C-16 


Small  diffraction  2  feet  below  the  ground  surface. 


6PR  Anomily  Number 

General 


LINBF 


Description 


F-14 


GPR  Anomaly  Number 

General 


General 


Thin  surface  layer  (1  to  3  feet  thick)  and  homo¬ 
genous  subsurface  sediments. 

Diffraction  at  5  feet  below  ground  surface  and 
some  ringing  and  reflectors  at  -20  feet. 


SINK  HOLE  SURVEY:  LINES  1~S 

Description 

Series  of  ns  lines  centered  on  sink  hole.  Lines  1 
and  2  show  generally  flatlylng  surface  layer  with 
homogenous  subsurface  stratigraphy.  On  lines  3. 4 
and  3  a  subsurface  reflector  can  be  seen  dipping  to 
the  southtowards  the  dam.  The  first  deflection  or 
rollover  occurs  at  the  sink  hole. 

SINK  HOLE  SURVEY:  LINES  6-10 

Series  of  ew  lines  centered  on  the  sink  hole  and 
oriented  at  right  angles  to  line  1-3.  On  lines  6.7 
and  8  a  pronounced  subsurface  escarpment  occurs 
approximately  30  feet  from  the  start  of  the  line. 
This  is  nearly  centered  on  the  sink  hole  and  the  ns 
line  number  3  (Lines  6.3  and  7.3  divide  these 
lines).  East  of  this  hinge  line  high  amplitude 
subsurface  refleaors  can  be  mapped  to  a  depth  of 
-13  feet.  On  line  9  and  10  a  shallow  subsurface 
depression,  centered  on  line  3.  extends  entirely 
across  the  survey  area.  On  line  10.  deeper  re¬ 
flectors  (-10  to  -13  feet)  below  this  depression 
indicate  the  presence  of  a  narrow  basin.  On  line  9 
a  small  subsurface  diffraction  can  be  seen  between 
cross  lines  1  and  2  at  approximately  3  feet  below 
the  ground  surface. 
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SUIIMAKT  AND  CONCLUSIONS 

Ground  penetrating  radar  was  successfully  used  to  map  the  subsurface  soil 
characteristics  along  4  transects  and  in  the  vicinity  of  a  sink  hole  at  the  Mill 
Creek  Reservoir.  The  following  summarizes  the  results  of  this  survey. 

*  In  general  the  subsurface  soils  along  the  primary  survey  lines 
(Lines  C  D,  E,  &  P)  were  quite  homogenous  and  uniform.  There  were 
occassional  subsurface  anomalies  but  the  data  does  not  clearly  in¬ 
dicate  that  these  features  on  the  radargram  are  subsurface  voids  or 
pathways  for  water  migration. 

*  In  the  sinkhole  area  several  features  or  anomalies  on  the  radar - 
gram  are  suggestive  of  possible  subsurface  voids  or  migration  path¬ 
ways.  On  line  9.  the  small  hyperbolic  feature  between  cross  lines  1 
and  2  is  a  typical  signature  of  a  void  like  feature.  On  line  10  the 
deeper  subsurface  depression,  between  cross  lines  3  and  <4,  is  pos¬ 
sible  evidence  for  a  zone  of  subsurface  migration. 

*  There  is  no  surface  ezpress  of  the  steep  subsurface  escarpment  that 
is  in  alignment  with  the  sinkhole,  Although  this  alignment  may  be 
coincidental,  several  small  hyperbolic  reflectors  beneath  the  toe  of 
the  escarpment  and  the  increase  in  apparent  reflectivity  suggests 
that  this  is  possibly  an  area  of  groundwater  leakage  and  piping  of 
sediment. 

*  There  was  no  evidence  of  conglomerates  or  bedrock  on  the  radar - 
grams. 
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PLATE  B8  Smknole  line  3. 
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PLATE  B9  Sinkhote  line  4 
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PLATE  B11.  Sinkhole  line 
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PLATE  812.  Sinkhole  line  6.5. 
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PLATE  Bt3.  Sinkho'e  line 
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Sinfchoie  line  T.5. 


LINE  8 


.  Sintehole  line  8. 


CROSS  LINES 


PLATE  816  Sinkhole  line  9. 
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PLATE  B17.  Sinkhole  line  10 
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Figure  C9.  Pole-dipole  survey,  line  C,  420  ft.  right/left,  ‘raw'  data 
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Figure  C11.  Pole-dipofe  survey,  line  C,  560  ft,  right/left,  'raw'  data  plot. 
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Figure  C12.  PoIe-dipoIe  survey,  line  C.  630  ft,  right/left,  'raw'  data  plot 
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Figure  Cl3.  Pole-dipote  survey,  line  C,  700  ft,  left,  'raw'  data  plot. 
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Figure  Cl 4.  Pole-dipole  survey,  line  C.  0  ft.  right,  smoothed  data  plot. 
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Figure  Cl5-  PoJe-dipole  survey,  line  C,  50  ft.  right/left-  smoothed  data  plot. 
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Figure  C19.  Po5e-dipole  survey,  litie  C,  250  ft,  right/left,  smoothed  data  plot. 
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Figure  C20.  Pole-dipole  survey,  line  C.  300  ft,  right/left,  smoothed  data  plot. 
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pole  survey,  line  C,  630  ft,  right/left,  smoothed  data  plot. 
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Figure  C26.  Pole-dipole  survey,  line  C,  700  ft,  left,  smoothed  data  plot. 


Figure  D1 .  Gravity  Data,  Line 
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Figure  D2.  Gravity  Data,  Line 


Figure  D3.  Gravity  Data,  Line 
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Figure  D4.  Gravity  Data,  Line 
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Figure  D5.  Gravity  Data,  Line 
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Figure  D6.  Gravity  Data,  Line 
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Figure  D7.  Gravity  Data,  Line 


Ik 


s 

^8  I  8  s; 

ilSoBell 


"a  s 

8  id  S  8 


%  8  15  8^8 

“2=S:isi 


£  3889883 


§8 

ii: 


I. IBs 

-Si?-! 


1  § 


8  8  S 


lUl 
■882 
8  8  8 


ill 

i  itf  :5  5 


imiini 


8888)(!8IC88 

iiiiiiiii 


888888888  I 

S  !C  8  if  8  8  »  a  S 
S8883(»SSS 

888888888 

iimim 


K.  «0  » 


8 

3 


•“  'O  >0 

ll'gi 

agsl 


O  'O  *#  •* 


ill 


%  ic  g 

K  c  e 


Se 


5  51 
K  C 


K  «  0< 
<0  iO  ««  «> 


OJ  o 
iO  O'  Ok 


i  a 

ft 

a 

n 

ft 

fM 

8 

ft 

a 

i  i 

% 

1 

1 

1 

1 

§ 

I 

8 

CS 

i  8 

8 

8 

3 

8 

8 

8 

8 

8 

:  « 

R 

S 

K» 

S 

C 

•O 

ft 

;  ft 

a 

a 

«n 

W\ 

m 

w% 

0% 

S 

8  8  8  8 
im 

»•  A»  M  k# 


8  8  8  8  8 
iim 

^O  K  10  0> 


8  ;a  8  s 

tiii 

8  e  8  e 


8  S  «  8 

O  «#  <«  Ok 

8  8  8  8 
°  8  R  K 


id  S  ft 

in 

o*  fk.  ^ 
Kt  PO  iO 
•O  K 

8  8  8 
?  *  * 

8  8  8 
ft  ft  3 


8  8  8  8 
o  e  o  o 

8  8  ft  if 

e  o  o  o 

8  8  8  8 
3  ft  ft  S 


8  8  8 
o  o  o 


IB  8 
r-i  o 

§  8 
ft  8 

S“ 

S  8 

n  ° 

8  8 
*  ° 

8  8 
d  d 


ft  ft  ft  SI  2 

d  ^  tO 


8  8  8 
•*  pi  K? 

R  C  I 


8  8 

i  si 


8!3!8ftft8S98 


el 


8  S  ft  ft  ft  8 

O  K  f*J  -p  K  ft 

»  ^  fM 

8  8  3  dt  $  8 

*>■  ‘f  ft  ft  R  5 

8  8  8  8  8  8 

°  S  ft  ft’  S  S 


8  ft 

~  a 

8  8 
5  ft 

8  8 
IS  ft 


8  8  8  8  8  8 
o  o  o  o  d  o 


8  ft  ft  ft  K  a 
°  ft  ft  R  8  S 


8  8  8  8  8  8 

s'gs'lsg 


8  8 
d  d 

R  8 

4  ? 
8  8 
s'g 


DIO 


Figure  D8.  Gravity  Data,  Line 


=  8.  j 

JH  J| 

t  '  ^  '  i  8 
H 


^  o  w  *rf.  «“. 

ieiiiiei 


^  ^  S  s  S ' 
E  c  g  S  4  * 


oordoo^^ 


soiiBsl  § 

8888888S  1 

88888888  T 

yiiieiii  i 


liiUtiiH 

,!  8tie9i::SS8 

!  j "  « 2  a  & « " 

i  8  e  t;  8  s  8  8  8 

i  i 

i 88888888 

Si®R«te»Bl‘‘ 

S  I 

I  8  8  8  8  8  8  8  8 

eedoC8000 

1  8  9  8  8  !s  8 
a  j  “'f'r’gcifcfifc 
s  * 

„!  8  8  8  8  8  8,  8  8 

iliiNslll 


Uiiim 


ddd^oooo 


88r:8i!!!^8i 


88888888 


88888888 


««  in  'fi  K  i* 


88bs;S^8l!8 

iiiHilii 

gRfttt8^88 

“«SH«  v" 

8  t;  S  9  i£  ti  8  8 

88888888 

“Ricrfija* 

8  8  8  8  8  ^  8  ^ 
ddddood9 

8Sl^i!r:8!i8 

88888888 

Uiiim 


Dll 


Figure  D9.  Gravity  Data,  Line 
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Figure  DIO.  Gravity  Data,  Line 
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Figure  DH .  Gravity  Data,  Line 
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Figure  D12.  Gravity  Data,  Line 
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Figure  013-  Gravity  Data,  Line 


siinssi  mt  vivf  coii 


my  SIHeI 

ygyiyasi 


S88gg8i88S 

gggggggggg 


gisgga^isRniSg 

°  !8 1  i  Hi  it  B  ll  S  ^ 

»  t  t  I  I  I  <  > 

8s;e;i«gs&:e@g 

®CCS?|SSS® 

lit* 

gggggggggg 

gggggggggg, 

•  eeoeigBOOO 

gK)gg|l;i!St3  98 
*»T‘?*'?7‘('t'rl’f 

I  88{^ggiggl(^ 


mmm 

"  i 

1 

b  It  g  K,  ^  g,  13 

lissgli 

IjyHsy 

1 

< 

1 

8  Ql  1  S  b  l{  8 

iimmi 

1 

i  1 

g  8  g  8  g  8.  g 

®  2  *i  g  «1  S  ® 

\  * 

'  iiiiiii 

e 

i 

gS!!5555 

i 

s 

i.  i|  g  g  ^  g  8 

D16 


Figure  Dl4,  Gravity  Data,  Line 
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Figure  015-  Gravity  Data,  Lina 
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Figure  Dl6.  Gravity  Data,  Line  F 
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Figure  017.  Gravity  Data,  Line 
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Figure  Di8.  Gravity  Data,  Line 
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Figure  D19.  Gravity  Data,  Line 


